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INDICATION TO AUTHORS 


Authors of papers concerning any of Smarandache type functions 
are encouraged to submit manuscripts to the Editors: 


Prof. C. Dumitrescu & Prof. V. Seleacu 
Department of Mathematics 
University of Craiova, Romania. 


The submitted manuscripts may be in the format of remarks, 
conjectures, solved/unsolved or open new proposed problems, notes, 


articles, miscellaneous, etc. They must be original work and 
camera ready [typewritten/computerized, format: 8.5 x 11 inches (= 
21,6 x 28 cm)]. They are not returned, hence we advise the 


authors to keep a copy. 

The title of the paper should be writing with capital letters. 
The author's name has to apply in the middle of the line, near the 
title. References should be mentioned in the text by a number in 
square brackets and should be listed alphabetically. Current 
address followed by e-mail address should apply at the end of the 
paper, after the references. 

The paper should have at the beginning an up to a half-page 
abstract, followed by the key words. 

All manuscripts are subject to anonymous review by two or more 
independent sources. 


A BRIEF HISTORY OF THE "SMARANDACHE FUNCTION" 


by Dr. Constantin Dumitrescu 
Department of Mathematics 
University of Craiova, Romania 


This function is originated from the exiled Romanian professor 
Florentin Smarandache. It is defined as follows: 


For any non-null integers n, S(n) is the smallest integer such 
that (S(n))! is divisible by n. 


The importance of the notion is that it characterizes a prime 
number, 1.@e.: 
Let p > 4, then: p is prime if and only if S(p) =p. 
Another properties: 

If (a,b) = 1, then S(ab) = max { S(a), S{(b) }; 
and 

For any non-null integers, S(ab) S S(a) + S(b). 

{All three found and proved by the author in 1979 (see [3], 15, 12- 
Ls, 65) :<} 

Ifin>> 1, then S(n) and n have a proper common divisor. 
{Found and proved by student Proddnescu in 1993: as a lemma needed 
to solve the conjecture formulated by the author in 1979 that: 

the equation S(n) = S(n + 1) has no solutions 
(see [3], 37, and [30]}).} 
EtG: 


Also, an infinity of open/unsolved problems, involving this 
function, provoked mathematicians around the world to study it and 
its applications (computational mathematics, simulation, quantum 
theory, etc.). 

Thus, the unsolved question: 

. n 

| 1 

Calculate lim!1 + ); —— - log S({n)!, (see [3], 29) 
Nee k=2 S(k) 

made by the author in 1979, has been separately proved by J. 
Thompson from USA in 1992 (see [18], 1), by Nigel Backhouse from 
United Kingdom in 1993 (see [25]), and by Pal Grgnas from Norway in 
1993 (see [51]) that this limit is equal to o. 

The author wonderred if it’s posible to approach the function 
(see [3], 1979, 25-6), but Ian Parberry expressed that one can 
immediately find an algorithm that computes S (n) in 
O(nlogn/loglogn) time (see [38], 1993). 

Some unsolved (by now!) other problems stated by the author in 
1979 (see [3], 27-30): 

a) To find a general form of the continued fraction expansion 
of S(n)/n, for all n22. 

b) What is the smallest k such that for any integer n at 
least one of the numbers S(n), S(ntl), ..., S(nt+k-1) is a 


perfect square? 
c) To build the largest arithmetical progression a,, a,, ..., 
a, for which their images by the function are also an 
arithmetical progression. 
Etc. 


In 1975 Smarandache was a student at the University of Craiova, 
and he was attracted by the Number Theory. He created and 
published a lot of proposed problems of mathematics in various 
scientific journals. He liked to play with the numbers 
Thus, in 1980 his research paper "A Function in the Number Theory", 
based on a special representation of integers, was published (for 
the first time) in <Analele Universitadtii Timisoara>, Seria Stiinte 
Matematice, Vol. 18, pp. 79-88, 
and was reviewed in <Zentralblatt fur Mathematik>, 471.10004, 1982, 
by P. Kiss, and in the <Mathematical Reviews>, 83c:10008, 1983, by 
R. Meyer. 

In 1988 he escaped from the Ceausescu’s dictatorship, spent 
almost two years in a political refugee camp in Turkey (Istambul 
and Ankara), and finally emigrated to the United States. 

Articles, notes, quickies, comments, proposals related to the 
Smarandache Function were presented to international conferences 
within the Mathematical Association of America or the American 
Mathematical Society at the New Mexico State University (Las 


Cruces), New Mexico Tech. (Socorro), University of Arizona 
(Tucson), University. of San Antonio, University of Victoria 
(Canada) etc. or published in <Octogon> (Sacele), <Gazeta 
Matematica> (Bucharest), <The Mathematical Spectrum> (UK), 


<Elemente der Mathematik> (Switzerland), <The Fibonacci Quarterly> 
(USA) etc. 

In 1992 Dr. J. R. Sutton from United Kingdom designed a BASIC 
PROCedure to calculate S(n) for all powers of a prime number up to 
a maximum. (see [26]) 

Jim Duncan from United Kingdom computed up to S$(1499999), the 
first million taking 50 hours in Lattice C on an Atari 1040ST. 
(see [17]) 

Also, John McCarthy from United Kingdom estimated that his 
machine would take several years to just calculate and store S (n) 
to disk for the entire range of n it can handle (0<n<2%*32), and 
using the compression detailed in ne1d9207.c at least 12 Gigabytes 
of disk space would be needed. It took about 3 hours for his 
program to work out that 3,303,302 pages (!) would be needed to 
list the full range of n and S(n). (see [15]) 

In 1993 Henry Ibstedt from Sweden used a dtk-computer with 
486/33MHz processor in Borland’s Turbo Basic and calculated $(n) 
for n upto 10° which took 2 hours and 50 minutes! (see [52]) 

A group of professors (V. Seleacu, C. Dumitrescu, L. Tutescu, I. 
Patrascu, M. Mocanu) and scientific students from the University of 
Craiova, having a weekly meeting, are doing research on the 
function and its applicability. 
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The Smarandache Function, S(n) 
{originated by Florentin Smarandache 
— Smarandache Function Journal, vol 
1, no 1, December 1990. ISSN 1053- 
4792) is defined for all non-null inte- 
gers, n. to be the smallest integer such 
that (S(n}}! is divisible by n. 

Note N! denotes the factorial func- 
tion, N!=1x2x3x...xN: for all positive 
integer N. In addition 0! = 1 by defini- 
tion. 

S(n} is an even function. Thatis, S{n) 
= S{-n) since if (S(n))! is divisible by n 
it is also divisible by -n. 

S(p} = p when p is a prime number, 
since no factorial less than p! has a 
factor p in this case where p is prime. 

The values of S{n) in Fig 1 are easily 
verified. For example, S(14) = 7 be- 
cause 7 is the smallest number such 
that 7! is divisible by 14. 

Problem (i} Design and implement an 
algorithm to generate and store/tabu- 
late S{n) as a function of n. 

Hint It may be advantageous to con- 

sider the STANDARD FORM ofn, vizn 
=ept: ps.......p4;-wheree=21,),,p,...-B, 
denote the distinct prime factors of n 
and a,,a,,....a, are their respective mul- 
tiplicities. 
Problem (ii) Investigate those sets of 
consecutive integers (i.i+1.i+2,...1+x} for 
which S generates a monotonic increas- 
ing (or indeed monotonic decreasing) 
sequence. 

Note For (1,2.3,4,5} S generates the 
monotonic increasing sequence 
0,2,3,4,5; herei=1&x=4. 

Lf possible estimate the largest value 
of x. 

Problem (iii) Investigate the existence 

of integers m.n,p,q & k with nem and 
q for which: 

either (A): S(m) + S(m+1) +....S(m+p) = 

S(n) + S(m+1) +....S(n+q) or (B): 

S(m} + S(m+1)}? +....S(m+p}? 

S(aF + S(n+1p +....9(n+p)? =k 


Problem (iv) Find the smallest integer k 


Republished from <Personal Computer World>, No.112, 


for which itis true that forall n less than 
some given n, at least one of: 

S(n), S(n+1)....S(n+k-1) is: 

A) a perfect square 

B) a divisor of k* 

C) a factorial of a positive integer. 

Conjecture what happens to k as n, 
tends to infinity: i.e. becomes larger 
and larger. 
Problem (v) Construct prime numbers 
of the form S(n) S(n+1)...S(n+k): where 
abcdetg denotes the integer formed by 
the concatenation of a.b.c.d,e,f& g. For 


example, trivially S(2) (3) = 23 which 


is prime, but _no_so_trivially 
5(14)S(15)S(16)S(17) = 75617, also 
prime! 


Definition An A-SEQUENCE is an 

integer sequencea,,a,,... with 1<a,<a,<... 
such that no a, is the sum of distinct 
members of the sequence (other than 
a.). 
Problem (vi) Investigate the construc- 
tion of A-CSEQUENCES a.,a,,.... such 
that the associated sequences 
S(a,),S{a,)....are also A-SEQUENCES. 

Definition The k® order forward fi- 
nite differences of the Smarandache 
function are defined thus: 

D,(x) = *modulus(S(x+1) - S(x)), 

D® (x) = D(D(...k-timesD,{x)...)} 
Problem (vii) Investigate the conjec- 
ture that D™(1)} = 1 or O for all k greater 
than or equal to 2. 

c.f. Gilbreath’s conjecture on prime 
numbers, discussed in ‘Numbers Count’ 
PCW Dec 1983. * Here modulus is 
taken to mean the absolute value of 
{ABS.}, modulus (y} = y if y is positive 
and modulus (y) = -y if y is negative. 

The following selection of 
Diophantine Equations (i.e. solutions 
are sought in integer values of x) are 
taken from the Smarandache Journal 
and make up: 

Problem (viii) [f m & n are given inte- 
gers, solve each of: 
a) S(x) = S(x+1), conjectured to have no 


420, July 


1992 (with the author permission), because some of the following 
research papers are referring to these open problems. 
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Republished from <Personal Computer World>, No.117, 


1992 (with the author permission). 


412, December 


solution 
b) S(mx+n) =x 
c) S(mx+n) = m+nx 
d)} S(mx+n) = x! 
e) S{x™) = x 
f) S(x}™ = S(x"} 
g) S(x) + y = x + S{y}. x&y not prime 
h) S{x) + Sv) = S {x+y} 
} 


i} S(x+v) = S{x}S(v 
j) SOx) = SONS) 


Review, July 1992 

The Smarandache Function: 

a first visit? 

This topic is certain to be revisited in 
the near future, and the lack of space 
available here will certainly be rem- 
edied on that occasion. Suffice it to 
report that Jim Duncan computed up to 
S(1499999}. the first million taking 50 
hours in Lattice C on an Atari 1040ST. 
In Problem (ii), no evidence fora largest 
value of x was found, while in Problem 
(vii) the conjecture was verified for the 
first 32.000 values of S(n). The very 
worthy prizewinner is John McCarthy 
of 17 Mount Street, Mansfield, Notts 
NG19 7AT, who has extensively inves- 
tigated the computation of S{n) up to 
2*; arriving at conclusions such as: 
‘several years of computing’, ‘at least 
12Gb of disk space’ and ‘3,303,302 pages 
ofoutput’. John’s concluding comment, 
‘Am I mad?’, is clearly answered NO! 
by examining his specimen pages of 
output including those relating to 10- 
digit values of n. Listings supplied. 
Details from Jobn directly upon request. 


ALGORITHM IN LATTICE C TO GENERATE S(n) 


by Jim Duncan 
9 Ryeground Lane 
Formby 
Liverpool 
L37 7EG ENGLAND 


Computer: Atari 1040ST 
Run time: ca. 50 hrs to generate S(n) for 1 000 000 numbers 


#incelude <stdio.h> 
#include <math.h? 


unsigned long int pst,ast,s,t; 


main () 
{ 


} 


long int n;3 
FILE *fp; 


printf ("input n (1 <2 n < 2 147 483 647)\n"); 

scant ("Zld"“, &n)3 

fp = fopen("PRNs*, “w"); 

fprintflfp, “n = Zid Stn) = Zld\n", nm, smaran(n))s 
fclose(fp); 


smar an (m) 
unsigned long int a; 


€ 


unsigned long int mst,p,a,fact; 


double rj 
if (m == 1) 
return(0); 
else ¢ /%* STANDARD FORM of m #/ 
r=mst = m3 
p = 1; 
fact = 13 
while ( ++#p <= sqrt(r)) € 
a = Q; 


while (mst ~% p == 0) ¢ 
mst = mst/p; 
att; 
+ 

r= mst; 

if (a >0o) ¢ 


pst = p; 
ast = a; 
£=s = 0; 
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3 


tors3(); /* find smallest factorial (t) with #/ 
/* p*a divsor ¥*/ 
if (t > fact) 
fact = ¢; 

} 
} 
if (mst > fact) 

fact = mst; 
return (fact); 
3 


tors! () 


€ 
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unsigned long int i; 


7/* test number is pst“ast #/ 

/* s is the difference between a factorial number t and ast#{pst-1) #/ 

/* s forms a pattern which determines the smallest value of t for which */ 
/* the test number is a divisor */ 


i = Q; 
while (++i < pst#pst && t-s < ast*¥(pst-1)) ¢ 
if (1 % pst == 0) 


Ss = s-pst+2; 
else 
st+5 
t += pst; 
} 
} 
tors2() 
a4 
unsigned long int i; 
tors1i);3 
i = 0; 
while (++i < pst*¥pst && t-s < ast#(pst-1)) ¢ 
if (i % pst == 0) 
Ss = s-3¥pst+43 
else 
S$ = s~-Zepst+3s 
t += pst; 
tors1()3 
} 
} 
tors3() 
€ 
unsigned long int i; 
tors2(); 
i = 0; 
while (++i < pst#¥pst && t-s < astx*(pst-1)) ¢ 
if (i Z% pst == 0) 
S$ = s-S*pstt+d; 
else 
Ss = s-4#pstt+5; 
€ += pst; 
tors2()5 
} 
3 


MONOTONIC INCREASING AND DECREASING SEQUENCES OF S(n) 


by Jim Duncan 


Froblem (11) 


Monotonic increasing and monotonic decreasing sequences of S(n) 
were investigated for x >= 6. 


First number (1) number of sequences S 

in range x = & x = 7 x = 8 x =9 
inc dec inc dec ;: inc dec inc dec 

1 - 499 999 735 83 7 10 0 2 1) ie) 

500 000 sare 999 999 80 7& 14 18 1 3 1 1 

1 Q00 000 - 1 499 999 7S 6&3 8 10 1 2 1 1 


There appears to be no evidence for a largest value for x. The sequences 
for x = 9 are shown in Results Table 1. 


The existence of sequences with the same first order finite differences 
was then considered eg: 


i = 440 Si) = 115 i= 5073 S(id) = 89 
i+i = 441 Siti) = 14; itl = 5074 Siti) = 59 
it2 = 442 Sti+2) = 17; it+2 = 5075 S(it+t2) = 29 


Apart from the initial quartet 2,3,4,5 all such sequences with 

i < 1 000 O00 are triplets. If the first order finite differences 
are multiples of 6 then the S(n) values appear to be prime numbers. 
The values are shown in Results Table 2. 
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RESULTS 


Sequences 


x 9333353593339 x 339339333393 x 9335393359393 98 x 


3939393329093 


inet enun ws 
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586951 
586952 
S84953 
586954 
586955 
586956 
586957 
586958 
586959 
386960 


9 


721970 
721971 
721972 
721973 
721974 
721975 
721978 
721977 
721978 
721979 


9 


1091150 
1091151 
1091152 
1091153 
1091154 
1091155 
1091156 
1091157 
1091158 
1091159 


9 


1473257 
1473258 
1473259 
1473260 
1473261 
1473262 
1473263 
1473264 
1473265 
1473266 


TABLE 


of S(n) 


S(n) 
S(n) 
S(n) 
S(n) 
Sin) 
Sin) 
Stn) 
Stn) 
S(n) 
S tn) 


| | | | 


S(n) 
S(n) 
S(n) 
Stn) 
S(n) 
S(n) 
Stn) 
Stn) 
S(n) 
S(n) 


Sin) 
S(n) 
S(n) 
S(n) 
S(n)> 
S(n) 
Si {n) 
Sn) 
S(n) 
S(n) 


S(n) 
Stn) 
S(n) 
Stn) 
Si (n) 
S(n) 
S(n) 
S(n) 
S(n) 
S(n) 


Htnt at www 


veut bt egntan in 


1 


586951 
73369 
21739 
9447 
1319 
1193 
1181 
1091 
677 
29 


73 

827 
907 
4067 
10939 
28879 
90247 
240659 
360989 
721979 


157 

709 
14351 
1607 
6271 
16787 
24799 
363719 
545579 
1091159 


1473257 
B467 
6323 
3877 
3533 
2239 
1999 
787 

557 

463 


RESULTS TABLE 2 
same difference = 1 nm = 4 S{n-2) = 2 S(n-1) = 3 Stn) = 4 
same difference = 1 nm = 5 S(n=-2) = 3 Stn-1) = 4 Stn) = 5 
same difference = 11 n = 18 S(n-2) = 6 S(n-1) = 17 Sin) = & 
same difference = 3S mn = 442 S(n-2) = 11 S(n-1) = 14 Sin) = 17 
same difference = 30 n = SO75 S(n-2) = 89 S(n-1) = 59 Sin) = 29 
same difference = 60 n = 6499 Sin-2) = 149 S(n-1) = 89 Sind = 29 
same difference = 48 n = $479 Sin-2>) = 127 S(n-1) = 75 Sins = Ti 
same difference = 36 nm = 8177 S(n-2) = 109 S(n-1) = 73 S(n) = 37 
same difference = 84 nm = 15717 Sin-2) = 214 S(n-1) = 127 S (ps ~ 43 
same difference = 168 n = 20468 S(n-2) = 379 Stn-1) = 721 Stn) = 43 
same difference = 210 nm = 22591 S(n-2) = 461 S(n-1) = 251 S(n) = 41 
same difference = 120 n = 35145 S(n-2) = Bit S(n-1) = 191 Sin) = 71 
same difference = 180 n = 59719 S(n-2) = 449 S(n-1) = 269 Sin) = 89 
same difference = 150 n = 67771 S(n-2) = 401 Stn-1) = 251 Stn) = 101 
same difference = 264 n = 73425 S(n-2) = 617 S(n-1) = ZS3 S(n) = 89 
same difference = 8640 n = 74005 S(n-2) = 1721 S(n-1) = @a1 Stn) = 41 
same difference = 24 nm = 82297 S(n-2) = 151 S(n-1) = 127 Sin) = 103 
same difference = 60 nm = 104669 S(n-2) = 251 S(n-1) = 191 Stn) = 131 
same difference = S30 nm = 111507 S(n-2) = 769 S(n-1) = 439 S(n) = 109 
same difference = 36 nm = 114427 S(n-2) = 199 S(tn-1) = 163 Stn) = 127 
same difference = 252 n = 120523 S(n-2) = 631 S(n-1) = 379 S(n) = 127 
same difference = 120 nm = 129978 S(n-2) = 389 Sin-1) = 269 Stn) = 149 
same difference = 952 n = 146004 S(n—-2) = 1973 Sin—-1) = 1021 S(n) = 69 
same difference = 600 n = 153520 Si(n-2) = 1301 S(n-1) = 7Ol Stn) = 101 
same difference = 12 n = 180482 S(n-2) = 47 S(n—-1) = S59 S(n) = 71 
same difference = 660 n = 181485 S(n-2) = 1429 S(n-1) = 769 Stn) = 109 
same difference = 60 nm = 189954 S(n-2) = S53 Stn-1) = 113 S (nm) = 173 
same difference = 90 nn = 192067 S(n-2) = 259 S(n-1) = 269 S(n) = 179 
same difference = 324 nm = 198697 S(n-2) = Ali S(n-1) = 487 S(n) = 163 
same difference = 336 n = 209752 S(n-2) = 939 S(n-1) = 5035 S(n) = 1467 
same difference = 228 nm = 227099 S(n-2) = 647 Sin-1) = 419 Stn) = 191 
same difference = 150 n = 231039 S(n-2) = 499 S(n-1) = 349 Stn) = 199 
same difference = 264 n = 253725 S(n-2) = 727 S(n-1) = 463 Stn) = 199 
same difference = 210 n= 266915 S(n-2) = 631 Stn—-1) = 421 Stn) = 211 
same difference = 648 n = 297638 S(n-2) = 1459 S(n-1) = 811 Stn) = 1463 
same difference = 2808 n = 306128 S(n-2) = 5669 S(n-1) = 28461 S(n) = 53 
same difference = 1320 nm = 324384 S(n-2) = 2749 S(n~1i) = 1429 Stn) = 109 
same difference = 18 nm = 326163 S(n-2) = 199 S(n-1) = 181 Stn) = 163 
same difference = 240 nm = 342965 S(n-2) = 719 S(n-1) = 479 Stn) = 239 
same difference = 34 nm = 3446390 S{(n-2) = 139 S(n-1) = 10S S(n) = 67 
same difference = 300 n= 386906 S(n-2) = 47 S(n-1) = 347 Sin) = 647 
same difference = 840 n = 409422 Sin-2) = 1861 S(n-i) = 1021 S(n) = 181 
same difference = 270 n= 440375 S(n-2) = B11 S(n—-1) = 341 Stn) = 271 
Same difference = 934 nm = 443450 S(n-2) = 2053 S(n-1) = 1117 Stn) = 181 
same difference = 120 n = 443850 S(n-2) = S09 S(n-1) = 389 Sin) = 269 
same difference = 792 n = 443969 S(n-2) = 1783 Sin-1) = 991 S(n) = 199 
same difference = 450 nm = 450043 S(n—-2) = 1151 Sin-1) = 7O1 Stn) = 251 
same difference = 306 n = 451215 S(n-2) = 883 S(n-1) = 377 Stn) = 271 
same difference = 210 nm = 460559 S(n-2) = 7O1 Si(n—1) = 491 Stn) = 281 
same difference = 240 nm = 464212 Si(n-2) = 761 Sin-1) = S21 S(n) = 281 
same difference = 360 nm = 470727 S(n-2) = 991 S(n-1) = 631 S(n) = 271 
same difference = 624 n = 473922 S(n-2) = 1481 S(n-1) = 8357 S(n) = 235 
same difference = 90 nm = 481779 S(n-2) = 449 Sin-1) = S59 S(n) = 269 
Same difference = 126 nm = 511688 S(n-2) = 131 S(n-1) = 257 S(n) = 383 
same difference = 672 n = 512894 S(n-2) = 1583 S(n-i) = 911 Stn) = 239 
same difference = 480 nn = 521946 S(n-2) = 1231 S(n-1) = 7351 Stn) = 271 
same difference = 714 nm = 531775 Stn-2) = 1667 S(n-!) = 9535 Stn) = 229 
same difference = 726 mn = 543455 S'(n-2) = 1697 S(n-t1) = 947 Sin) = 241 
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same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 
same 


difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 
difference 


nudtrrkt vu uau nt tun na 


io vot on 


506 
5352 
1176 
3444 
390 
840 
336 
612 
78 
264 
498 
S40 
1356 
102060 
S214 
1590 
252 
552 
2850 
540 
240 
120 
1122 
264 
684 
744 
310 


n 

n 
n 
n 


93390395 
il 
Hou non 


303 


393939 


ia] 
fy 
nr 
n 
n 


a 


ia 


9335939 


565187 Sl(n-2) 
= 3574498 S(n~2) 
= 586272 S(n-2) 
592537 S(n-2) 
609871 S(n-2) 
629508 S(n-2) 
681077 S(n-2) 
705793 S(n-2) 
724319 S(n-2) 
726827 S(n-2) 
731179 S(n-2) 
731746 Sin-2) 
778837 S(n-2) 
792475 S(n-2) 
803427 S(n-2) 
810451 S(n-2) 
= 837969 S(n-2) 
= 898783 S(n-2) 
930311 S(n-2) 
941057 S(n-2) 
943553 S(n-2) 
97353546 S(n-2) 
997443 Sin-2) 
1026550 S(n-2) 
1028985 S(n-2) 
1042162 S(n-2) 
1053175 S(n-2) 


“oudouwou 


= 919 S(n-1) = 613 S(n) = 307 

= 1381 S(n-1) = 829 S(n) = 277 
= 2549 S(n-1) = 1373 Sin) = 197 
= 6971 S(n-1) = 3527 Stn) = B3 

= 1091 S(n-1) = 7O1 S{n) = 311 

= 1931 S(n-1) = 1091 S(n) = 251 

= 1009 S(n-1) = 673 Sin) = 327 

= 1531 S(n-1) = 919 Stn) = 307 
467 S(n-1) = 389 Stn) = Sit 

= 881 S(n-1) = 617 Sin) = 253 

= 1327 S(n-1) = 829 S(n) = 331 

= 829 Sin-1) = $99 Gin) = 359 

= 2939 S(n-1) = 1583 S(n) = 227 
= 2311 S(n-1) = 1291 Stn) = 271 
= 4391 S(n-1) = 2377 Stn) = 162 
= 3391 S(n-1) = 1801 S(n) = 211 

= 883 S(n-1) = 631 S(n) = 379 

= 1471 S(n-1) = 919 S(n) = 367 
= S851 S(n-1) = 3001 Stn) = 151 
= 1459 S(n-1) = 919 S{(n) = 379 

= BB1 S(n-1) = 641 Stn) = 401 

= 619 S(n-1) = 499 S(n) = 379 
= 2551 S(n-1i) = 1429 S(n) = 307 
= 947 S(n-1) = 683 Stn) = 419 
= 1747 S(n-1) = 1063 Stn) = 379 
= 1861 S(n—-1) = 1117 S(n) = 373 
= 1429 S(n-1) = 919 Stn) = 409 


ON THE CONJECTURE D,™)(1)=1 or 0 for k>=2 


by Jim Duncan 


Problem(vi1) 


For the first 32 000 Sin)‘s the conjecture that De **? (1) = 1 ofr O for 
k >= 2 is true. The ratio of the number of ones to the number of zeros 
appears to be approximately { for large values of k. The results are 
shown in Results Table 3, 
The true differences S(x-1) - S(x) were calculated and the kt order 
differences D.* **? (1) were found to increase rapidly with increasing k. 
For large values of k ( > 100) the ratio Dy "<*? (1) /Do *Sk-272 (1) is 
approximately equal to -2. Some values are shown in Results Table 4. 


RESULTS TABLE 3 
Di=2. 
Ratio: Di000 = O sti/stO = 1.112051 
k oat ee) total 1’s D2000 = 0 sti/stoO = 1.056584 
total 0O’s D3000 = 0 sti/stO = 1.075423 
D4000 = O sti/stO = 1.015625 
DS0O0O0 = O sti/stO = 1.010661 
D6000 = 1 sti/stO = 0.991039 
D700Q = 9 sti/stO = 0.990048 
D8000 = O sti/stO = 0.977014 
D9000 = 0 sti/stO = 0.987412 
D190000 = oO sti/std = 0.998201 
Di1000 = i sti/stO = 1.011154 
Di2000 = 1 sti/stO = 1.010556 
D13S000 = 1) stl/stO = 1.015036 
D14000 = 0 sti/stO = 1.018601 
D1S000 = 1) sti/stO = 1.012748 
D16000 = O sti/stO = 1.904134 
D17000 = 0 sti/stO = 1.005308 
D18000 = O sti/std = 1.004789 
D17000 = 1) sti/stO = 1.003903 
D20000 = 1) sti/stoO = 1.004711 
D21000 = 0 sti/stdO = 1.008129 
D22000 = 0 sti/stO = 1.004830 
D23000 = 0 sti/stO = 1.004620 
D24000 = 0 sti1/stO = 1.003590 
D25000 = 1 sti/stO = 1.004571 
D26000 = 1 sti/stO = 1.001001 
D27000 = 0 st1l/stO = 1.001260 
D28000 = 1 sti/stO = 1.004080 
D29000 = 0 sti/stO = 1.006018 
D50000 = 1 sti/stO = 1.005415 
DS1000 = 0 sti/stO = 1.006408 
D32000 = 1 sti/stO = 1.004699 
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RESULTS 


k Do) 


2.2444216+288 
—4.496445+288 
9.0099146E+288 
-1.805740E+289 
3.619671E+289 
-7. 257009E+289 
1.455178E+290 
-2.9183466E+290 
5. 853595E+290 
-1.174247E+291 
2.3558275+291 
-4. 72681 92+291 
9. 484823E+291 
-1,903346E+292 
=. 8194685E+292 
~7. 645702E+292 
1. 5384526+293 
~3. OB7569E+293 
6.1964545+295 
-1. 24357154294 
2. 495458E+204 
-5.007457E+294 
1.004734E+295 
~2.015791E£+295 
4.043841E+295 
-8.111306E+295 
1.4267846+296 
~3.262162E+296 
6. 540525E+296 
-1.3111305+297 
2.627848E+297 
-5. 2458856+297 
1. 055006E+298 
-2.113228E+298 
4.231969E+298 
-8.473041£+298 
1. 696031E+299 
~3.394086E+299 
6.79053 1E+299 
-1.358233=+300 
2.71601 3E+300 
~5. 4296288E+300 
1.085180E+361 
-2.1682757E+301 
4.3311296+301 
-B8. 64911984301 
1.726721E+302 
-3.446291E+302 


= 6.876396E+302 
~1.371668E+ 3023 





TABLE 4 
Ds (1) 

DY (1) 
ratio = -2.002974 
ratio = -—2,003387 
ratio = ~2.003787 
ratio = -2.004170 
ratio = -2.004535 
ratio = -2.004881 
ratio = -2.005204 
ratio = -2.005504 
ratio = -2.005778 
ratio = -2,.006026 
ratio = -2.006246 
ratio = ~2.006437 
ratio = -2.006598 
ratio = -2.006728 
ratio = -2,.006827 
ratio = —-2,.904894 
ratio = -2.506929 
ratio = -2,006932 
ratis = -2.a06904 
ratio = -2.006643 
ratiq = -2.o06751 
ratio = -2,.0066729 
ratio = -2.006476 
ratio = -2.006293 
ratio = -2.006081 
ratio = -2.005842 
ratio = -2.005576 
ratio = -2,005283 
ratio = -2.004966 
ratio = -2,004625 
ratio = -2.004262 
ratio = -2.003877 
ratio = -2.003473 
ratio = -2.003049 
ratio = -2,002608 
ratio = -2.002151 
ratio = -2.001679 
ratio = -2.001193 
ratio = -2.000695 
ratio = -2.9000184 
ratio = ~1.999667 
ratio = ~1.999139 
ratio = ~1.998604 
ratio = -1.998063 
ratio = ~1.997516 
ratio = -1.996966 
ratio = -1.996415 
ratio = -1.995858 
ratio = -1.9952303 
ratio = ~-1.994748 


A Simple Algorithm to Calculate S(n) 
by John C. McCarthy 


Introduction 


This short paper first outlines an "“obvious" algorithm for 
calculating S(n) (the smallest integer m such that m! is divisible by 
n). Doubtless, there exist more subtle and efficient algorithms. I 


hope some readers will devise these and enlighten me concerning them 
through this journal. 


This is followed by a small scale investigation of the efficiency 
of the algorithm. 


Then there is a short discussion of a simple way of reducing the 
space required for storage of all S(n) for ranges of n. The storage 
space required for S(n) for all n which my routines can handle is 
considered. 


Heavily commented listings of an implementation of the algorithm 
in "C”, sample output and timing data are included to help illustrate 
the algorithm. 


he Algorith 


The algorithm is described in detail at the start of the header 
file "“S(n).H". Together with "S(n).c", this forms all the code 
necessary to implement the algorithm. Note that, for the S(n) function 
to work correctly, the function make_primes() must first be called from 
the main program. 


The code for printing S(n) and timing the routines has been 
omitted. These activities are both implementation specific and easily 
done. They are therefore left as an exercise for the interested 
reader. 


The algorithm hinges on finding the prime factors of on. 
Improvements on how this is done will most benefit its efficiency. 


To be practical, the given implementation of the algorithm only 
works for 0<n<232. However, the algorithm is generally applicable to 
any non-null integer. 


Tables of S(n), constructed using the routines of "S(n).C", for 
the largest 2000 permitted mn are included.’ My paging 
routines are rather elaborate. Using them (without printing!), it took 
2.4 hours to discover that 3,745,708 pages, as tightly packed as those 
shown, would be required to print S(n) for all O<n<23?. 


Efficiency of the Algorithm 


In a letter to R. Muller (about computing the Smarandache 
Function, July 19, 1993), Ian Parberry (editor of <SIGACT News>, 


1 
For the smallest 4800 numbers, sea Ibstedt’s table (pp. 43-50) of 19 
this current journal. 


Denton, Texas) expressed that one can immediately find an algorithm 


that computes S(n) in O(nlogn/loglogn) time ('A Brief History of the 
"Smarandache Function"' by Dr. Constantin Dumitrescu, Department of 
Mathematics, University of Craiova, Romania). Disappointingly, a 


little analysis of the accompanying timing data on my TI85 advanced 
scientific calculator reveals that my algorithm is somewhat worse than 
this. 


Trying to fit the version 2 timing data to various O(f(n)), I 
obtained the following results (x=3355443200 and 10(0(x+99)-O(x-100)) 
is calculated for comparison with the last entry of the version 1 
timing data): 


10(O0(x+99)-O(x-100) 


O(f£(n)) | Correlation | O(232-1) 
(milliseconds) 


Coefficient (years) 














8909 









O(n) 0.9928879 0.6092 
O(nlogn/loglogn)| 0.9944006 0.7906 11827 
O(n/n) 0.9997756 24.2 469178 


O(nvyn) fits the version 2 timing data best, although the time it 
predicts for the last entry of the version 1 timing data is almost 3 
times too large. Hence, I assume the time complexity of my algorithm 
is a little better than O(n/n). 


As a rough upper limit on the time my program (on my 20MHz 368DX 
PC) would take to calculate S(n) for all 0O<n<2°?, let us assume that 
every value of n requires as much time as each n in the range of the 
last entry of the version 1 timing data (= 159111/199/10 = 79.9553 ms). 
In this “worst case”, it would take 10.882 years. O(nvyn) time 
complexity predicts more than twice this value, which is a measure of 
how pessimistic it is. 


I would welcome a more rigorous analysis of the time complexity of 
my algorithm as I presently lack the necessary expertise. 


Simple Compression of Stored S(n) 


Without compression, each S(n) would be stored as a 32-bit (= 4 
bytes) value. Hence 2°* bytes (= 16 Gigabytes) would be required to 
store S(n) for all O<n<232. 


This requirement can be reduced considerably if we use the high 
bit of each each byte of each value to indicate if it is the last byte 
of the value. If the bit is set it means that further byte(s) are 
required and if it is reset it means that the byte is the last byte of 
the current value. This means that only 7 bits of each byte are used 
to form the numerical part of the value. Assuming that, as with Intel 
format, the values are stored low-'byte' {actually 7 bits) first, here 
are some examples: 


i) 127 requires seven bits and so just one byte (with high bit 


reset to indicate no further bytes). 
ii) 16,000 requires 14 bits. So it is stored as two bytes. The 
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first is 0 (16,000 mod 128) + 128 (to set the high bit 
indicating there is more to come). The second is 125 (16,000 
div 128) (with high bit reset to indicate no further bytes). 
This reads simply as 0 (with more to follow) + 128*125 (no 
more to follow). 

iii) A number stored as the three bytes 57+128, 93+128 and 125+0 
would similarly represent: 
57 + 93*128 + 125*128*128 = 2,059,961. 


The largest numbers that can be represented by a given number of 
bytes is thus as follows: 


127. 


1 byte can code up to 27-1 = 

2 bytes can code up to 2?4%-1 = 16,383. 

3 bytes can code up to 27%-1 = 2,097,151. 

4 bytes can code up to 27%-1 = 268,435,455. 

5 bytes can code up to 2°5-1 = 34,359,738,355 (or 8 times the 


largest unsigned long). 


For small values of n, the savings are considerable (400%). 
However, even large n often have small S(n). 


Using this technique to compress all S(n) calculated for some 
ranges of n (each range was also stored), I obtained the following 
results: 


range of n compression time taken size size after pkzip 
(seconds) 
1 without 4.5 40,008 19,836 
-~10,000 with 4.7 15,749 15,267 
2,147,478,648 without 827.3 40,008 33,729 
-2,147,488,647 with 842.4 33,541 30,836 
4,294,957,296 without 1,066.2 40,008 34,320 
-4,294,967,295 with 1,085.1 34,330 31,634 


The results indicate that this compression is a little better than 
pkzip's (a commercial file compression utility). Application of pkzip 
to a pre-compressed file also gives a slight improvement. 


Assuming that the savings shown for the middle range of 10,000 n 
are the average of all ranges of 10,000 n, using my compression 
together with that of pkzip would permit storage of S(n) for all 
O<n<22? in about 3.0836*2°? = 12.3344 Gigabytes. So look out for sets 
of 19 CD-ROMs with all your favourite numbers on them! 


2lst November 1993 
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/* 


af 
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S(n).H 


(c).1993.11.13.John.C.McCarthy 
"S(n).h" 


Example Implementation of A Simple Algorithm to Calculate S(n), 
The Smarandache Function: 


Because there are more people familiar with C than with C++, this module 
has been written entirely in C (apart from "//" style comments). The 
module was compiled using Borland C++ version 3.1. 


For efficiency, n is constrained to the limits of an unsigned long. Hence, 
0 <=n <= 2732 - 1 (= 4,294,967,295). (""" represents exponentiation). 
Although catering for n of vast magnitude is possible, it imposes heavy 
storage and processing overheads. The range of an unsigned long therefore 
seems a reasonable compromise. 


The algorithm depends on the most elementary properties of S(n): 
1) Calculate the STANDARD FORM (SF) of n: 


In SF: n = +/-(pli7al)*(p27a2)*...*(pr°“ar) where pl, p2,...pr denote the 
distinct prime factors of n and al, a2,...ar are their respective 
multiplicities. 

2) S(n) = max(€ S(pl7al),...,S(prar)]. 

3) S(p7a), where p is prime, is given by: 
3.1) at=p == S(p°a) = p*a. 
3.2) a>p == S(p7a) = x < p*a. In this case, fortunately rare, x 


is the smallest integer such that p appears as a factor in the 

list of all integers > 1 and <= x at least a times. Let the no. 

of times p appears as a factor in the list of all integers > 1 

and <= y be fly, p). Then: 

fly, p) = Elintly/(p7i))] for i>0 while yo=(p7i). 

Hence, x is the smallest integer such that f(x, p)>=a. 

Note that between succesive integer multiples of p there are no 

integers which have p as a factor. The trick here is to look for 

the largest multiple of p (call it ¢c), such that f(p*c, p)<=a 

(so that x = p*c, if f{p*c, p) = a, else x = p*(ctl1)): 

3.2.1) ¢ = a-2 (largest possibilty for c since f(p*(a-1), p)>=a 
when a>p (Note: f(p*(a-1), p)=a is not sought for slight 
performance gain)). 


3.2.2) 2 = f£(p*c, p). 

3.2.3) While(z>a): 
3.2.3.1) d = no. of times p appears as a factor of p%*c 

= (no. of times p appears as a factor of c) + 1. 

3.2.3.2) c = c-l (next largest possibility for c). 
3.2.3.3) 2 = 2-d (= £(p%c, p)). 

3.2.4) If(z<a), x = p*(ctl1). 

3.2.5) Else X = p*c. 


To calculate the prime factors of all 32-bit n requires use only of primes 
< (2716) (i.e. all primes expressible as an unsigned short integer). This 
is because any factor of n remaining after division of n by all its prime 
factors < (2716) is simply a prime. Since there are only 6542 16-bit 
primes, the program first creates a list of these (which only takes about 
4 seconds on my 20 MHz 386DX PC) so that they never have to be 
recalculated, thus saving much time. 


S(n).H 


#define PRIMES16 6542 // The number of 16-bit primes 

#define MAX_SFK 9 /* max. distinct primes in the SF of n. The smallest 
number with more than 9 distinct primes is the product of the 10 smallest 
primes (= 6,469,693,230), which is substantially more than the largest 
integer expressible as an unsigned long. Hence, 9 distinct primes are 
more than ample. 

/ 

typedef unsigned long u_long; 

typedef unsigned int u_int; 

typedef enum {false, true} boolean; 


struct SF_struct { 


int sfk; // no. of distinct primes 
u_long sfp(MAX_SFK]; // the distinct primes 
int sfalMAX_SFK]; // vespective multiplicities 


}; 


extern u_int prime[PRIMES16+1]; // list of all 16-bit primes 
// plus terminating zero. 


void make_primes(void); // construct list of all 16-bit primes (prime[]). 
_ // Must be called before calls to getSF() or S(). 
void getSF(u_long n, struct SF_struct *SF); // calc. SF of n and store in SF 
u_long S(u_long n); // cale. S(n) 
u_long Spa(u_long p, int a); // calc. S(p”a) where p is prime 
int f(int x, int p); /* the number of times the prime p appears as a factor 
in the integers from 1 to x inclusive. This function is only called from 
Spa(p, a) when a>p with x=p*(a-2) (refer to item (3) of algorithm outline 
above). Max value of (a) occurs when p is a minimum, n is a maximum and 
(p*a)=n. So, (27max(a))=max(n)=(2732)-1. Hence max(a)<32. So, x<60 
when (a) is at its max. Max value of p (and x) occurs when a=ptl and 
(p°a)=max(n). So, max(p)7(max(p)+1)=(2732)-1. The upshot is that 
max(p)=9 when a=10. Hence, max(x)=72. This explains why it is safe for 
x, p and the return value of f(x,p) to be passed as ints. 

x/ 
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/* 


x/ 


S(n).c 


(c).1993.11.13.Jonn.C.McCarthy 
"S(n).c" 


Example Implementation of A Simple Algorithm to Calculate S(n), 
The Smarandache Function: 


This is the code for the module. Refer to "S(n).h" for details. 


#include "S(n).h" 


u_int prime[PRIMES16+1]; // allocate storage for list of all 16-bit primes 


// plus terminating zero. 


void make_primes(void) 


{ 
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u_int *pp; // ptr to last prime so far of prime list 
u_int *tp; // ptr to current test prime 


u_int p; // number being tested for primality 

pp=prime; // point to start of prime list 

*pp=2; // set first prime to 2 

*++pp=3; // set second prime to 3 

p=5; // next possible prime. N.B. p is kept odd so that trial 


// division by 2 is unnecessary. 


whibe(true) { // infinite loop!: 


tp=prime+l; // point to first odd test prime 
// whilst test prime <= Vp: 
while(((long) *tp)*(*tp)<=p) { 
if(!(ps*tp)) { // If current test prime divides (is factor of) p: 


pt=2; // try next odd number 
if(p<*pp) { // done when p overflows: 
*++pp=0; // terminate list 
return; 
tp=primet+l; _// point to first odd test prime 
else ++tp; // Else point to next test prime 
} 
// no prime <= vp divides p so p must be prime: 
*++DD=p; // so store it next in the list 
pt=2; // try next odd number 
if(p<*pp) { // done when p overflows: 
*++pp=0; // terminate list 
return; 


void getSF(u_long n, struct 
{ 


S(n).c 


SF_struct *SF) 


u_int *pp; // ptr to current prime 
u_long r; // ‘vesidue' of n remaining for factoring 


SF->sfk=0; // no. of distinct prime factors discovered 


r=n; 


pp=prime; // point to start of prime list 


// whilst current prime <= vr and prime list not exhausted: 
while(((long) *pp)*(*pp)<=r && *pp) { 


if(!(r*pp)) ¢{ 
SF->sfp(SF->sfk]=*pp; 
SF->sfa(SF->sfk]=1; 
r/=*pp; 
while(!(r%*pp)) { 
SF->sfa(SF->sfk] ++; 
r/=*pp; 


SF->sfk++; 
} 
++pp; 
} 


if(ro1) f{ /f 
SF->sfp(SF->sfk]l=r; 
SF->sfa(SF->sfk]=1; 
SF->sfk++; 

} 

} 


// if current prime is a factor of r: 
// store current prime as next prime of SF 
// set its multiplicity to 1 
// ‘divide out' current prime 
// while current prime factors r: 
// increment multiplicity 
// ‘divide out' current prime 


// increment count of distinct prime factors 
// next prime 
If n contains prime > 2716: 
// store it as last prime of SF 


// set its multiplicity to 1 
// increment count of distinct prime factors 


25 


uw. 


{ 


} 


S(n).c 


long S(u_long n) 


struct SF_struct SF; // to store SF of n 


int sfi; // index of current term of SF of n 

u_long Sn; // current guess at S(n) 

u_long x; // S(current term of SF of n) where it might exceed 
// current value of Sn. 

if(n==1) return 0; // special case 

getSF(n, &SF); // cale. and store SF of n 


// First guess at S(n) is S(p a), where p is the largest prime in the SF 

// of n and a is its multiplicity. This pre-empts the calculation of S(p*a) 
// for the remaining terms where, as is likely, p*a for these terms is <= 

// this initial guess (since S(p°a) <= p*a always): 

sfi=SF.sfk-1l; 

Sn=Spa(SF.sfp[(sfiJ],SF.sfa(sfil]); 


while(sfi>0) { // while more term(s): 


sfi--; // next term 
if(SF.sfp[(sfil]*SF.sfalsfiJ>Sn) { // if this term may have larger S(pa): 
x=Spa(SF.sfp(sfij,sF.sfa({sfil); // cale. it 
if(x>Sn)} Sn=x; // if new max., update Sn with it 


dy 


} 
return Sn; // That's all folks! 


u_long Spa(u_long p, int a) 
{ 


} 
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// Refer to item 3) of the algorithm description in S(n).h. 

int c; // largest multiple of p such that f(p*c, p) <= a (eventually!) 
int 2; // f£(p*c, p) 

int m; // used to calc. no. of times p appears as factor of c 


if(a<=p) return p*a; 


c=a-2; 
2z=f£(p*c, p); 
while(z>a) { 
// ain items 3.2.3.1) and 3.2.3.3) of algorithm description is implicit 


// nere: 
aes 
M=C~--; 
while(!(m&%p)) { // while p divides m: 
Z--; 
m/=p; // ‘divide out' factor of p fromm 


} 


if(z2<a) return p*(c+1); 
else return p*c; 


S(n).c 


int f(int x, int p) 


{ 
int k=0; // count of appearance of prime p as a factor in the integers 


// from 1 to x. 
int xdp; // successive divisions of x by p 


xdp=x/p; 

while(xdp>0) { 
k+=xdp; 
xdp/=p; 


return k; 
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SMARANDACHE FUNCTION, 


Scad, 


\S of a plus: 


AN e 
4294965296 268435331 
4294965506 S531 


4294965316 19737413529 
42949653286 2147482663 


4294965336 21517 
4294965346 3919 
4294965356 14969 
4294965366 aig 
4294965376 $527 


4294965386 2147482693 
4294965396 357913783 





4294965406 235289 
4294965416 3179 
4294965425 816223 
4294965436 63983 
4294965446 178793 
4294965456 8329 
4294965466 10275037 
4294965476 157 
4294965486 102261083 
4294965496 Sogsi 
4294965506 17903 
4294965516 357913793 
4296965526 21474682763 
4294965536 5835551 
4294965546 199798 
4294965556 4948117 
4294965566 13174741 
42794965576 16268809 
4294965586 19701677 
4794965596 1073741399 
4294965606 1543973 
4294965616 315231 
4294965626 682513 
42734965636 9177277 
429496564 2937733 
4294965656 20731 
4294965666 906211909 


4294965676 10737414619 


4294965686 195225713 
4£294965696° 2729 
42394965706 1959439 
4294965716 220889 
4294965726 217166193 
4294965736 536870717 


4294965746 2147482873 


4294965756 2818219 
4294965766 10276031 
4294965776 268435361 
4294965786 4286393 
427949657596 a337 
4794965806 23831 
4294965816 761 
4294965826 420829 
4294965836 359231 
4294965846 1018247 
4294965856 7129 
4294965866 107869 
4294965878 63113 
42949658856 2147482943 
4294965896 19289 
4234965306 6871 
4294965916 12064511 
4294965926 z256049 
4294965936 89478457 
4294965948 7451137 
4294965956 5417 
4794965966 1191061 
4294365976 76695822 
4294965986 7430737 
4294965996 6907 
4294966006 36398017 
4294966916 15102 
4294966026 Z79511 
4294966036 $6512711 
4294966046 238001 
4294966056 4519 


4294966066 2147483033 


42943966076 103613 
82394966086 79536409 
4294966096 20543 
4294966106 721474383053 
4294966116 131441 
4294966126 2749 
4294966136 536870767 
4294966146 19933949 
4294966156 4858559 
4294966166 15673059 
AZ94966176 481067 
4294966186 126613 
4294966196 1012009 
4294966206 1063637 
4294966216 14510021 
42794966226 16007 
4294966236 1121987 
4294966246 2147483123 
4294966236 38347913 
4294966266 163543 
42946966276 859681 
4294966286 165191011 
4294966296 04867 
42794966306 251197 
4294966316 26188819 
4294966326 1571 
42949663356 17159 
42794966346 195225743 
4294966356 56893 
4294966365 283571 
4294966376 267249 
4294966386 3s5t 
4294966396 32783 
4294966406 116098 


i 

163837 
4294965307 
458397461 
64s101 
11392481 
4294965347 
158597 
30971 
830587 
21367987 
4060419 
2347 
625997 
9965117 
1530091 
65573 
4294965457 
14364433 
1194041 
4294965487 
14387 
72091 
99882919 
252645031 
26881 
107077 
20355287 
1431655189 
1618909 
1398 

218941 
118362 
4294965617 
1828423 
682933 
209623 
26372 
27111617 
Z28azsz73 
1431655229 
Savar 
394353 
1286303 
849312 
42394965737 
887573 
4294965757 
4294965767 
1404961 
244157 
$5823 
41357 
6771138 
2129 

7883 
4294965847 
237173 
1186127 
36749 
4294639658587 
841321 
364691 
TSZ217 

31is 
42394965937 
z63 
38693387 
4294965967 
4294965977 
SGE7OS9I12% 
5956957 
4294966007 
14737 
35495587 
2235797 
163003 
2213 
613566531 
477218453 
42794966087 
44278001 
669937 
3271109 
73583 

$21 

340519 
186737659 
153801: 
4294966177 
4234966187 
47279 

2113 
42794966217 
830887 
4294966237 
S65693 
477218473 
3$5573 
87652373 
35391 
8294966297 
2458481 
5940479 
116080171 
42394966337 
68174069 
3761 
4294966367 
1307 
186737669 
21013681 
1431655469 


z 3 

2799847 60492469 
1194373 101377 
102261079 273564667 
7255009 26723 
$320S21 22199183 
e911 90407 

31909 108907 
4925419 1431655123 
134573 279857 
3567247 30034723 
35461 27012362 
1354807 329393 
221338997 4oi13 
12109 6040739 
79536397 7473 
36097 364607 
367067 423191 
7615187 1822217 
2147482739 1935541 
14639 973253 
718827583 390451409 
29483 3804221 
4699087 12236369 
37193 4294965529 
88919 3296213 
1073741387 23339 
109383937 12569 
67108837 397057 
3081037 1491655193 
zsa3 23549 
9617603 22605082: 
336870701 923063 
2957 89611 
46634409 2459 
2147482819 14312655213 
89478451 119747 
147827 4234965659 
1139557 4603593 
FIS8Z7G13 4294965679 
400949 20338427 
17S3)0 477228421 

21131 75431 
2147482859 8ss4i1i 
601873 8783161 

179 7517 
63161261 650851 
7877 204522179 

zsit 24793 
1162637 603989 
84407 1811 
$7649 919 
24317 1160860157 
349241 2165893 
3579138719 4799 
S3S5319 4294965839 
36599 14316552853 
72211 148192271 
1673741467 33353 
SOG78S277 477218431 
2179 2081903 
2147482949 613566557 
1073741477) 110127331 
187759 1848% 
1574491 SZ182 
1622 130150483 
3933129 4294965948 
6221 138547289 

2591 8699 
88389 1321 
26188817 13258 
33332 42389 
178956917 60492479 
4153133 2579 
258173 254879 
65075243 47197429 
1342137689 273617 
21467483023 539069 
3224449 9695183 
182129 91382257 
3701 1431655363 

45137 4294966999 
1973741527 35083 
2147483059 33294311 
257863 4294966129 
2147483069 252645067 
8873897 3697 
77647 390451463 
905347 251 
74051141 2237 
188397571 80737 
2008871 7548272 
Si162219 $7847 
43481 237409 

3737 446323 

95287 1199071 
462421 28676809 
171491 12744707 
1073741567 381673 
7WLS8Z27713 7363 
163779 16417 
38839333 332712 
357913859 127663 
108191 6074917 
536879791 159072827 
3411889 325697 
77813 39403361 
2147483179 73350287 
3832029 401867 
113025431 399451489 
151637 2464123 
17459213 399941 
6047 41698703 


for n=#4294965296 to n=4294967235 


4 
24316551 
5174657 
16025399 
37427 
12632251 
83641 
2339 
135137 
1199712 
7535027 
1952257 
$167 
z3a6osis 
321239 
13421767 
286339103 
70019 
49699 
18233 
18673763 
1227133 
235197 
33687069 
324509 
2468371 
2371683 
28057 
20452217 
113907541 
429496559 
19884) 
9283 
3147 
251699 
2551 
g5899313 
17459213 
429496567 
101875 
47771841 
42949657 
61356653 
1556147 
79099 
16319039 
1Z1843 
1677722 
7229 
4933 
3538 
21474823 
143165527 
57327 
33669 
1626881 
721843 
214748293 
407879 
107376147 
2161 
14316553 
429496591 
103643 
181913 
2214748297 
4651 
627919 
991969 
690503 
229 
13% 
13159 
$939 
2225371 
120581 
9544369 
27889399 
10475527 
1033 
4013987 
azs 
233549 
221712 
$74249 
1213 
acoag 
53687077 
13015048 
93396883 
119737 
1022611 
1633067 
10853 
199673 
246271 
33607 
13463 
356661 
761287 
7863 
261412 
919633 
5052 
430789 
1614649 
28633109 
107053 
429496637 
24253 
273043 
1342177 
3497 


61 
9 


$s 

3226871 
66569 
4919777 
4294965331 
40139863 
1431655117 
4294965361 
LIB974I1 
204522161 
354493 
3778513 
144611863 
135107 
28793 
1331089 
117321 
4294965461 
88873 
22031 
2754949 
835759 
2017 
832519 
103391 
e529 
627549 
2293s 

13709 
4294965581 
9739151 
4294965601 
35527 
7495577 
9137 
42794965641 
62245875 
16582879 
4294965671 
18013 
4294965691 
3373 

132867 
4294965721 
494413 
139150477 
226050873 
187553809 
477218419 
12377423 
41698697 
70524659 
38008547 
4294965821 
@51641 
4294965841 
954649 
11096103 
613566553 
186737647 
366997 
252645053 
4294965911 
405683 
97553 
asacs 
79259 
30634727 
4294965971 
21431655327 
61333 
42949669001 
243769 
49031 
116080163 
329927 
3256229 
54366659 
3853 
12321767 
11392483 
85733 
4554577 
4294966121 
1759 
1328477 
BE41783 
S168431 
B2139 
3904514671 
21367991 
3521233 
2245147 
1748053 
4294966231 
24239 
1431655417 
$501 

96263 

4073 

Sil 
72796039 
477218479 
226050859 
19489 
2706343 
135347 
246681 
1431655457 
908603 
33818633 
3557 

22067 


8 
49941457 
967 
21487482651 
3741259 
2286989 
1228537 
2147482681 
SS7913782 
39107 
268435337 
70867 
672349 
17747791 
S772803 
308857 
221313 
415937 
169681 
718943 
7253 
16909313 
193583 
715827587 
32789 
7306641 
89478449 
7561 
51829 
193748173 
8801159 
2147482801 
307 
2147482811 
170327 
784039 
100547 
1187357 
5396523901 
505409 
130261 
94399 
72727 
118273 
23057 
1689601 
29173 
89669 
1000691 
836573 
7433 
2102559 
19627 
297383 
$36870729 
2147482971 
32749 
399391 
268435367 
280312 
6869 
2147482951 
87767 
25799 
B31713 
639703 
3947579 
306783283 
13597 

git 
117709 
715827667 
1073741503 
197018679 
1439 
33368827 
3343 
24683713 
949223 
303617 
S3S7913841 
3067832393 
7064089 
14035837 
242873 
S$8o040083 
4153 
288679 
16067 
55063669 
24403217 
6317 
$235371 
302466241 
$36870779 
23743 
1109 
16393001 


97612873 
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? 8 

231547 337867 
4294965313 113025403 
22967729 2117833 
75577 853889 
5895613 379691 
2545919 715827559 
23719 118423 
104755253 641 
42934965383 6628033 
S749619 2147682697 
26349481 Siise9 
4294965413 34301 
204522163 39727 
746561 3003 
2617 28871 

13997 2465537 
138547273 536870683 
4001 5549051 

84869 1073741371 
30509 28661 
637519 294337 
7033 306783251 
142799 1073741381 
6381821 12132671 
211831861 64319 
179623 17599 
77933 219304599 
20507 2986801 
12781 2949839 
72796027 2970239 
332943907 40459 
42794965613 78487 
330381971 4409 
1302689 2147482817 
7027 4549 
6400843 160033 
28071671 20929 
4294965673 333685819 
4294965683 1360889 
45613 113025413 
13054607 291907 
2238127 183% 
249287 153391633 
4294965733 2147482857 
7309 6701 
142837 2147482877 
B2LSS 1073741441 
19792469 65075239 
@2Z2Z3) 5536870723 
4294965793 306783271 
556559 14512 
1431655271 17471 
81037091 20719 
2395829 2386909213 
48661 1737351 
346117 4303573 
4564061 235565273 
65173 2147482937 
70309 1091 
620923 8949783 
9633 1909 

7829 7388279 
125941 351931 
zz43 438276183 
1290167 536370743 
23718359 1494421 
21431655321 241453 
311749 1390253 
74aa7l 31867 
5282861 S276371 
1322749 38053 
471301 21318283 
6537239 S29981 
3637 353691 
4294966043 3579133837 
238649 1263225312 
351263 28627 
4294966073 11987 
46182431 7110871 
1081583 8893 
340573 2657 
167269 69273647 
ZE3IS33 1823 
665783 31122943 
792283 4793489 
4294966153 2147483077 
4294966163 S78213 
6661 9629969 
817933 97987 
14965039 191 
28313 2862629 

6719 21503 

1361 3263439 
1739557 4029049 
4294966243 1073741561 
2116080163 763957 
29671 821 
2217379 2147483137 
35537 58169 
G721387 74051143 
35747 1254371 
3904651483 15230377 
49367429 1600211 
613566619 49941469 
232645079 178956931 
$4133 200381 
1762399 1367 
4294966373 170557 
159072829 268435339 
1057613 126322541 
210900L 357913867 
149333 1679033 


9 
2798023 
286331021 
246342659 
8589393067 
4621 
1493 
33629 
1307 
858993077 
86357 
977239 
24683 
10105801 
1453457 
858993089 
2207 
168331 
4079 
8863 
137593 
72679 
18841? 
1847297 
122713301 
106747 
8919397 
78090283 
1931 
337 
50529007 
3335361 
21859 
3061 
3133283 
40153 
40961 
838993133 
3Ot1 
12347 
7880671 
286331047 
838993143 
29153 
859853 
21553333 
6270023 
4457 
65497 
16943 
286331053 
8753 
122713309 
19088737 
49993 
39409 
733 
7933 
15953 
3843491 
780907289 
8S8993181 
3443987 
172798637 
3886847 
19927 
122713313 
858993193 
457 
SO513 
8589931399 
32497 
1340083 
19219 
11603 
1128769 
277363 
16231 
171798643 
68791 
13634813 
1433 
780902393 
11453243 
2340581 
27709459 
286331077 
1227133139 
171798647 
86331 
365063 
199765387 
Z8G6331081 
15618059 
812671 
15070057 
9323 
29620457 
179 
858993257 
373475353 
40904441 
313547 
197243 
397 
8389937269 
838993271 
25733 
701229 
73727 
19991 
838393281 
858993283 


SMARANDACHE FUNCTION. S(n), for n=4294965296 


\ 


Ay 
4294966416 
4294966426 
4294966436 
4294966446 
22794966456 
4294966466 
4294966476 
4294966486 
4294966496 
42946966506 
4294966516 
4294966526 
4294966536 
4294966546 
4294966556 
42942966566 
4294966576 
4294966586 
4294966596 
4294966606 
4294966616 
4294966626 
4294966636 
4294966646 
4294966656 
4294366666 
4294986676 
4294966636 
42794966696 
4294968706 
4294966716 
4294966726 
4294966736 
42794966746 
4294966756 
4294966766 
4294966776 
4294966786 
4294966796 
4294966806 
4294966816. 
4294966826 
42794966856 


2334388838 


4294966866 
&Z94966876 
4294966886 
4294966896 
4294966996 
4294966916 
4294966926 
47949663936 
4Z94966946 
42794966956 
4234966966 
4294966976 
4294966966 
4294966996 
4294967006 
4294967016 
4294967026 
4294967036 
42949676046 
4294967056 
4294967066 
4294967076 
429439670856 
4294967096 
42794967106 
4294967116 
4294967125 
4294967136 
4294967146 
4294967156 
4294967166 
4294967176 
4294967186 
4294967196 
4296967706 
4294967216 
4294967226 
42949672356 
4294967746 
4294967256 
42949677266 
4294967276 
42949672386 


S of n plus: 


6 1 2 
1877 813566631 48677 
1564081 4294966427 7017919 
373571 17674739 so9g 
3912627 4294966447 1172207 
26361 288427 6567227 
306783319 130150499 2848121 
38891 6294966477 3308911 
69273653 613566641 178956937 
7895159 49367431 2147483249 
234467 416623 1153 
306872 49891 1951 
2231 477218503 16777213 
BS21759 482093 2147483269 
265291621 252645091 16193 
733883 1129957 8910719 
34781 65563 463219 
1405422 793 37633 
33040089 7TS35OZ9L 133169 
383 aaa 9727119 
15661 38835159 29826157 
870131 105323 764501 
79536419 45751 12064513 
879671 330382049 7839 
21474839323 1281697 1080223 
17449 4294966657 4121849 
S56487 423934966667 152239 
757 251771 49941473 
17959 390451517 2383 
336870837 16582883 21692751 
21474833553 975907 4492643 
219304631 3607 919 
S2377643 64103981 13765919 
1020667 30460757 11864549 
115821 180461 33653 
1073741689 Ssiises 1299143 
321047 191 27751 
1844917 9133 99223 
1490273 148102303 34019 
1789 433835903 21474833599 
761 99682949 32173 
3257 1516049 359171 
29417581 i3z29 9502139 
: 48623 rerererer 13449521 
14 sé 
7387403 15722088 121988 


191449 2999 1073741717 
YO73741719 4294966877 5076793 
306783349 333213 23633 
1489) «6330382069 1246363 

3639 178267 357913909 

7307 16928669 82339 
7FiSS27SZL 4294966927 372357 
18121 7403451 42107519 
20069939 $9390191 §$6$12723 
13813907 353693 23087 
70067? 226950893 59652319 
67108859 163487 2147483489 
12558383 369851 6072887 
119159 4294966997 26107 
11483869 65519 5835553 
38149 795217 14412641 
23917 218541 32537629 
2445881 1413283 3355967 
16111 39153 41319 

S347 252645121 83459 

15761 132967 zs5998S59 
119304641 23761 $239571 
214748635439 €294967687 1132639 
76695841 1431655699 2147483549 
1383 2135737 2609 

83841 10058471 754297 
2147483563 2184 85639 
13147 138547327 1467863 
93368ES51 I3sggo0izt 2251631 
1073741789 1932039 2147483579 
1262483 3229 1s0% 
B7481 104755297 317159 
252293 204522247 3449 
Z2$21 4294967197 1296907 
1809169 91387281 18539 


3780781 1431655739 $2377649 
ZZO037 252645131 106311907 


417961 825473 29389 
ZB1ii 65543 LST7A7A 
$9652323 $2613 2147483629 
367531 37349 295387 
361651 1861711 218046081 
715827681 65539 2089 


to n=4294967295 


3 

8088449 
$5569 
15913 
1137137 
£72787 
17970571 
7001 
3455323 
4999961 
477218501 
1241933 
613566647 
16086017 
3023 
330382043 
$11123 
138547309 
12739 
6388627 
22253713 
4294966619 
23469763 
42794966839 
284831 
139361 
214309 
18979 
36709117 
103961 
169761353 
70164163 
208889 
32292983 
14316555383 
198391 
4294966769 
477218531 
851399 
17821438 
294522279 
1766749 
4294966829 
83609 


3481333 


2909 
911689 
1931 
154423 
429429669083 
641327 
8263 
186737693 
47197439 
397351 
97289 
zsa791 
1588963 
16330673 
1597 
1599053 
42949670238 
zos4s 
477218561 
17355197 
2408843 
LLO1Z7361L 
613566727 
96553 
1431655703 
15053 
161471 
159072857 
2944249 
3156023 
179383 
62929 
42794967189 
461382443 
4s8saag 
11959 
69779 
704039293 
153966421 
asgi7 

35809 
4294967279 
1039 


4 

251 

2939 
761251 
365162 
202021 
143165549 
681 
429496643 
36767 
81083 
8259551 
443237 
3963947 
12271333 
64373 
201547 
214748329 
4772185) 
32587 
1523 

6531 
1297573 
5822 
2202547 
88627 
429496667 
6053 
21401 
a7119 
163165557 
431 
144563 
713582779 
17179867 
47389 
477218533 
16519103 
7459 
3579139 
429496681 
214748341 
35081 
349753 


2366384 


293 
268439543 
3491843 
6135667 
2029 
3254561 
5438667 
19522577 
$3923 
364397 
762% 

347 

93081 
4294967 
143165567 
180007 
11608019 
ra 53 

soe9 
214748353 
1431655869 
1O7374177 
429496709 
362 
2252177 
236507 
463319 
214748357 
858939343 
647 
184571 
34981 
143165573 
31033 
25264513 
25411 
23854733 
23327 
#111 
32771 
115861 
127 
27605091 


3 

30899039 
204822211 
4294966441 
2669339 
19611719 
12763 
330352037 
159072833 
813566643 
1868189 
345727 
2210873 
Z3I078S1 
1431655517 
119083 
235577 
21577 
4294966591 
465781 
119127349 
1693 
312839 
15671 
4294966651 
4294966661 
1237 
24826397 
109453 
73350293 
2280917 
2194 

84349 
40129 
41421723 
$9417 
116089183 
2387419 
362979 
93229 
1072937 
13899563 
913872773 
38351 


908333 


330332067 
1431855627 
8311 
113783 
1431655637 
seo0gz1se 
104755291 
8297 

9817 
447439 
1431655657 
4294966982 
3427747 
110127353 
22453 
365747 
$683 
10710641 
22582749 
7121 

6$551 
417839 
21467 
181967 
4294967111 
102871 
12898 

1163 

77249 
4294967161 
651839 
3583 

82791 
11576731 
477218579 
39403369 
4294967231 
14631655747 
1664071 
281 

3449 

g951 
4294967291 


6 

4804213 
377 
614093 
78301 
126322543 
6468323 
321143 
S$6512717 
sais 

443 
3984199 
4493 
T1S827757 
536870819 
$2377641 
5107 
22139003 
222953 
z20369 
5088823 
21262221 
719 
74051149 
1073741663 
20063 
8659297 
86477 
357913891 
38569 
336870639 
482689 
(397 
195225761 
14913079 
46411 
107021 
TISS2Z7797 
336870549 
93503 
32779 
26561 
349941 
20347 


1030383 
1345541 


38940033 
1073741723 
715827817 
70051 
BSS55711 
37943 
306783353 
536870859 
217643 
449 

202397 
323027 
4106087 
21913997 
5549053 
536870879 
96703 
357913921 
4007 
1Z2Z016%2 
192261121 
3203 
66377 
182423 
12859183 
10737461783 
2476913 
16381 
3agig 
357913931 
1935225781 
23342213 
238609289 
292493 
813749 
Bi91 
145681 
1$27371 
19867 
536870909 
795659 
333 


7 
390451499 
389848 
1431655482 
3979521 
998382941 
1949 
2176871 
6806603 
1431655501 
390487 
9565627 
32611 
613566649 
4294966553 
617359 
186163 
4294966583 
21367993 
208261 
7993 

92347 

8597 

14549 
799361 
736023 
1033189 
120307139 
54366667 
138547313 
1431655571 
693971 
8sgaas 
539459 
69463 
15176561 
24265349 
6938557 
330382061 
1431655601 
4294966813 
29078399 
477218537 
isz2Ssser 
782333 
35963 

1327 
138053 

git 

S979L 
13687 
453199 
12744703 
4294966943 
1431655651 
87652387 
390451543 
1431655661 
18341 
$54977 
375467 
17401 
3764213 
9851 

37Z02 
768193 
160s6ar9 
3904513553 
21332367 
68174081 
88463 
17321 
gseo2z1 
4294967143 
33818639 
48861397 
9739 

14767 
23333 
26029 
751687 

3089 

3557 
619139 
1431655751 
64103989 
395959 
1226783 
464773 
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8 

7213 

12541 
41539 

199 
2485513 
2147483237 
1073741621 
7866239 
536870813 
1236573 
10631 
2941753579 
149713 
238609253 
530767 


19SZ2S737 
63997 
TIS8ZT7779 
1073741671 
1028981 
12782639 
70% 
29259277 
4079 
36209 
79841 
736689 
306783341 
178481 
10683997 
a4ai 

4021 
25679 
27059 
273383 
*437233 
47731 
€11993 
27531839 
30647 
22133 
2797 

3272 
35204647 
915% 
2147483477 
1073741741 
VLS8Z7829 
78341 
2147483497 
1323569072 
12707003 
14128181 
3392543 
346368312 
2879 
$64533 
1779191 
342283 
21691753 
$197 
675097 
81031 
186689 
7321 
70853 
381707 
21474683587 
6197 
7134437 
97612891 
65981 
1726273 
318949 
45083 
862097 
282657 
6118187 
46684427 
2147483647 


3 
37266219 
636763 
656221 
$039 
8589393293 
133382 
2951867 
8589932399 
6151 
26407 
7469507 
40123 
56509 
32443 
16207421 
364753 
9439487 
179593 
286331107 
490573 
16901 
2219621 
4567 
15877 
26030101 
120223 
5867 
286331113 
218629 
2221 
4241 
9324643 
838993349 
35671 
27709463 
15618061 
134523 
7S311 
317323 
88129 
61291 
8589933967 
1373 
653557 
4517 
2290649 
2212139 
3651611 
35443799 
858993383 
36191 
21839 
858993389 
683 
393853 
12413 
41177 
133 
84358597 
858393403 
10141 
858993407 
11767033 
22284889 
19976591 
791693 
12917 
3953903 
9199 
2863391241 
14753 
129347 
376119 
315293 
858993433 
7283 
44371 
435210181 
286331147 
13209 
231223 
22025473 
69001 
4259 
95443717 
171798691 
122713351 
65537 


Timing of S(n).c Module for Calculation of Smarandache Function, version 1 





Time taken to calculate S(n) depends on how easy it is to factor n. Less 
time is required if n has "small" prime factors. So, in the following table, 
the values of n shown are the mid-points of ranges (n-99 thru nt+99). Times 
shown are for calculating S(n) for all integers in each range 10 times over: 


n time (ms) 





100 268 

200 308 

400 345 

800 387 

1600 432 

3200 490 

6400 571 
12800 661 
25600 766 
51200 919 
102400 2450 
204800 4036 
409600 5670 
819200 7977 
1638400 10423 
3276800 13004 
6553600 16302 
13107200 23438 
26214400 29642 
52428800 37011 
104857600 50330 
209715200 62363 
419430400 77888 
838860800 108179 
1677721600 158480 
3355443200 159111 


30 


Timing of S({n).c Module for Calculation of Smarandache Function, version 2 


"Time to n" is the time taken to calculate S(n) for all n <= that shown. 
"Time add.” is the time taken to calculate S(n) for all n> previous n and 
<= current n. All times are in milliseconds (as per version 1): 


Time add. 

































































50000 18223 18223 
100000 66763 48540 
150000 139191 72428 
200000 229634 90443 
250000 335252 105618 
300000 452539 117287 
350000 579419 126880 
400000 715146 135727 
450000 859963 144816 
500000 1012335 152372 
550000 Li71 221 158886 
600000 1336899 165678 
650000 1508825 171927 
700000 1686808 177983 
750000 1870023 183215 
800000 2058983 188961 
850000 2252457 193473 
990000 2450892 198435 
950000 2653620 202728 
1000000 2860734 207115 
1050000 3072049 211314 
1100000 3288502 216454 
1150000 3509106 220603 
1200000 3733965 224860 
1250000 3962171 228206 
1300000 4194158 231987 
1350000 4429331 235173 
1400000 4668560 239229 
1450000 4910513 241953 
1500000 5155601 245088 
1550000 5404652 249051 
1600000 5656512 251859 
1650000 5911306 254794 
1700000 6169686 258380 
1750000 6431383 261697 
1800000 6696172 264789 
1850000 6963206 267034 
1900000 7232974 269768 
1950000 7505412 272438 
2000000 7779763 274351 
2050000 8056579 276816 John McCarthy 
2100000 8336442 279863 17 Mount Street 
2150000 8620053 283611 Mansfield 
2200000 8905641 285588 Notts. 
2250000 9194727 289086 NG19 7AT 
2300000 9486449 291722 United Kingdom 
2350000 9780105 293655 
2400000 10076920 296815 
2450000 10375202 298282 
2500000 10676383 301181 31 





Mike Mudge pays a 
return visit to the Florentin 
Smarandache Function.? 


The originator of this function. 
Florentin Smarandache, an Eastern 
European mathematician. escaped from 
the country of his birth because the 
Communist authorities had prohibited 
the publication of his research papers 
and his participation in international 
congresses. After spending two years 
in a political refugee camp in Turkey, 
he emigrated to the United States. 
Robert Muller of The Number Theory 
Publishing Company, PO Box 42561, 
Phoenix, Arizona 8508, USA, decided 
to publish a selection of his papers, 
mmencing with The Smarandache 
Function Journal. Vol 1, No 1, Decem- 
ber 1990. ISSN 1053-4792. 
PCWreaders mav have met this func- 
tion before, in Numbers Count -112- 
July 1992, where a very encouraging 
response was generated. This article 
(February 1993] is complete in itself so 
don't worry if vou have filed the July 
issue! It mav be thought that those 
readers who attempted the previous 
problem-set will have an unfair advan- 
tage. However, it must be realised that 
no Numbers Count problems are com- 
pletely original so previous work within 
agiven subject area is alwavs a possibil- 
ity and the prize is awarded using ‘suit- 
able subjective criteria’ anyway, so 
please have a go and submit your re- 
sults. however trivial they may seem to 
yourseif. 
Definition For all non-null integers, n, 
the Smarandache Function, S(n), is 
defined to be the smallest integer such 
that (S(n))! (The Factorial Function with 
argument S(n),) is divisible by n. e.g. 
S(18) = 6 because 6! is divisible by 18 
but 1!....5! are not. 
Problem (0) Design and implement an 
algorithm to generate and store/tabu- 
late S(n) as a function of n upto a given 


NG. <7 
Hint It may be advantageous to con- 
sider the STANDARD FORM ofn, vizn 
= ep’, pi: p& ......p%s where e=+1, and 
P,.P,-P,---P, denote the distinct prime 
factors of n and a,, a,, a,...a, are their 
Tespective multiplicities. 

NOTE S(n) is an even function, by 
which is meant S{-n) = S(n). 

Problem (i) Using either graphical or 
finite difference technique (i.e. the con- 
struction of difference tables etc) or 
indeed anything else that comes to 
mind. address the following questions: 
(a) Is there a closed expression (for- 
mula) for S(n}? 

(b) Is there a good asymptotic expres- 
sion for S({n)? (By which is meant a 
formula, which although never (in gen- 
eral) exact, becomes a better and better 
approximation to S(n) as n becomes 
larger and larger.) 

Problem (ii) For a specified non-null 
integer m. under what conditions does 
Sin} divide the difference n - m? 
Problem {iii) Investigate the possible 
integer solutions, (x.y.z] of S(x*} + S(y*) 
= S(z*) for any no greater than or equal to 
1, e.g. examine the solution (5,7,2048) 
when a = 3. 

(It can be proved that an infinity of 

solutions exist for anv such n-value.) 
Compare with Fermat’s Theorem re. 
xe+y? = z8 
Problem (iv) Investigate the possibility 
of finding two integers n and k such 
that the LOGARITHM of S{n*) to the 
BASE S(k*) is an integer. 
Problem (v) Recall that ‘Gamma’ de- 
fined as the limit as a tends to infinity 
of (1+1/2+1/3 + 1/4....+1/n - log(n) } 
exists, is known as Euler’s Constant 
and is approximately 0.577. 

Investigate the possible existence of 

‘Samma’ defined as the limit as n tends 
to infinity of (1 + 1/S(2) + 1/S(3)....4+1/ 
S(n)} - log(S{n)}. 
Problem (vi) Find the number of PAR- 
TITIONS of n as the sum of S(m) for 
2<msn. See PCW August 1989 and 
February 1990 for other problems in- 
volving PARTITIONS of n. 
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Review of ‘Numbers Count -118- 
February 1993: a revisit to The 
Florentin Smarandache Function? 
This produced a number of ‘quite 
powerful’ responses. As a note of 
related interest, the latest publication 
of Fi.Smarandache is ‘A Numerical 
Function in Congruence Theorv’, 
Libertas Mathematica (American 
Romanian Academy of Arts and 
Science) vol 12, 1992, pp 181-185, 
Arlington, Texas. 

Pal Gronas of Norway submitted 
theoretical results on both problems 
O & {v). However, the clear winner 
this month is a former regular respon- 
dent, now retired. Henry Ibstedt, 
Glimminge 2036, 280 60 Broby, 
Sweden. Henry used a dtk-computer 

with 486/33MHz 
processor in Borland’s 
Turbo Basic. S(n) upto 
105 took 2hr 50min. 
He completed a great 
deal of work on all 
problems except {vi): details of numer- 
ical results and conclusions available 
from Henry or myself to interested 
readers. What about problem {vi}? 


A NOTE ON S(p‘) 


by Pal Grgenas 
Enges Gate 12 
7500 STagRDAL 
Norway 


t=1 


Problem (0). If []£, p{' is the prime factorization of n, then it is easy to verify that 


k 
S(n) = S(T] pi) = max{ 5(pi*) is 
t=1 
From this formula we see that it is essensial to determine S(p’), where p is a prime and 


r is a natural number. 
Legendres formula states that 


k - = 
ni= II»: Lim=iln/?i V 


t=1 
This formula gives us a lower and an upper bound for S(p"), namely 
(1) (p—1)r+1 < S(p’) S pr. 
It also implies that p divides S(p”), which means that 


S(p") = p(r —1t) for a particular 0 <2 < (=: 
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A proof of the non-existence of ”Samma’”. 


by Pal Grgnas 


Introduction: If []#., p' is the prime factorization of the natural number n > 2, then 
it is easy to verify that 


k 


S(n) = S([[ pi") = max{ S(p*) JE. 


t=1 


From this formula we see that it is essensial to determine S(p”), where p is a prime and 
ris a natural number. 
Legendres formula states that 


(1) = [Tp 2amr 


=1 


The definition of the Smarandache function tells us that S(p’) is the least natural num- 
ber such that p”|(S(p7))!. Combining this definition with (1), it is obvious that S(p") must 
satisfy the following two inequalities: 


Me 
| aa 
7 W 
ays, 
eae 


(2) So [set Seel<r < 


k=1 


> 
I 
ro 


This formula (2) gives us a lower and an upper bound for S(p”), namely 


(3) (p—1)r+1 < S(p’) < pr. 


It also implies that p divides S(p"), which means that 





rah). 


S(p") = p(r —1) for a particular 0 <i < [ - 


”"Samma”: Let T(n) = 1 — log(S(n)) + D2, ch Hy for n 2 2. I intend to prove that 
limyoo T'(n) = 00, i.e. "Samma” does not exist . 
First of all we define the sequence p; = 2, pp = 3, ps = 5 and p, = the nth prime. 
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Next we consider the natural number p™,. Now (3) gives us that 





S(pk) < pk Wie {l,...,m} and VRE {I,...,n} 
4 
it ee ste 
S(pk) ~ pk 
4 
m n l m nr l ™ 1 nr l 
Lesh 2 YLaet (Sp) (La) 
y 
(4) vay ~ >) 
aa (kK) Ne? SK 


since S(k) >0 for all k > 2, p} <p? whenever a<m and b<n and p? = p? if and 
only if @a=c and b=d. 

Futhermore S(pi,) < pm, which implies that —log S(p%,) > —log(pmn) because log z 
is a strictly increasing function in the intervall [2,00). By adding this last inequality and 
(4), we get 


T(p,) = 1—log(S oy +S a 08(Pm n)+(S-)- (04) 
= ka Pk k=1 
4} 
T(r) > cian (n = Pm) 
4 
Tipe) 2 142 (—logam+ 5) + (-24 E>) (FF) 
4 
dim, Tem) 2 1+ 2- im (—logmn+ 30g) + Jim, [(-2+ 0 —)-(¥9)] 
k=1 Dk k=l 
: . <1 
= 1+ 2- lim (- wae p + dal(- 2450): (a) 


II 


Pm 1 
1+ 2y + jim (- 2+ a —)- i pum, m (>> 2 (y = Euler’s constant) 
k=1 


oO 


It 


since both )4_, ¢ and Y4_, =, diverges as t — oo. In other words, lima. T(n) = 
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A BASIC PROCedure to calculate S(p%i) 


by John Sutton 
16A Overland Rd. 
Mumbles, SWANSEA 
SA3 4LP, U. K. 


Integer function of a single variable S%(N%) 
S is the least integer such that S$! 1s divisible by N. 


Obviously for a prime S(p)=p since this is the least factorial to 
include p. 


It is easy to see that for two primes pl>p2 S(pl*p2)=pl since this 
factorial is necessary to include pl and already includes p2. This 
generalizes to the product of any number of primes. 


In fact it generalizes to the product of relatively prime numbers nl 
and: n2. S(n1*n2)=Max(s(nl),S(n2)). 


Therefore we can simplify the general case to: 


s (Uni*pi) =Max (S(ni*pi)) 


All we need now is a way of calculating S for powers of primes. 


Start with the inverse problem: for a given factorial and a given 
prime what is the maximum power of the prime included ? 


Consider p=2. All even numbers contribute a factor, all multiples of 
4 contribute another, all multiples of 8 contribute yet another 
_..ete. So the answer is got by summing succesive DIV 2 results 


(DIV p in general). 


Returning to the calculation of S. To do this for a single N would 
require factorisation of N first. A program to calculate S_ for all 
integers up to N can avoid this by doing powers of 2, then powers of 
3 and their products with powers of 2 then powers of 5 etc. 
Calculating S for all powers of a prime up to a maximum is 
straightforward. A BASIC PROCedure is attached. The main program 
requires some care and I have not been able to finish in time. 


LOREM TEST PROC TOCALC S(P“I) 
20: 
30: 
40: 
50: 
6CINPUT"UP TO",N% 
TJODIM SPP%(100) 
80DIM NPP%(100) 
QOINPUT"WHICH PRIME",P% 
LOOPROCSpp(P%,N%) 
LIOFOR I%=0 TO 100 
L2OPRINT SPP%(1I%),NPP%(1I%) 
L30NEXT I% 
140GOTO 60 
LSOEND 
160DEF PROCSpp(P%,N%) 
1L70I%=1 
L8ONPP%(0)=1 
190SPP%(0)=1 
2007%=1 
210PI%=0 
220REPEAT 
230PI%=PJ%+P% 
240X%=FNinvSpp(P%,PJ%) 
250REPEAT 
260SPP%(I1%) =PJ% 
270NPP%(1%) =P%*NPP%(I%-1) 
2801 %=1%+1 
29O0UNTIL I%>X% 
Z00T%=I%+1 
31LOUNTIL NPP%(I%-1)>N% 
320ENDPROC 
330DEF FNinvSpp(P%,N%) 
340LOCAL S%,T% 
350S%=0 
360TS=NS 
370REPEAT 
380T%=T% DIV P% 
390S%=S%+T% 
400UNTIL T%<=1 
410=S% 


FOR VALUES UPTO N 
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Henry Ibstedt 


Glimminge 2036 
280 60 Broby 
Sweden 
The Florentin Smarandache Function S(n) 
PCW February 1993 
Problem (0) 


A program SMARAND has been designed to generate S(n) up to a preset limit N (N up to 
1000000 has been used in some applications). The program requires an input of prime 
number up to VN. Initially the program caculates S(p,*)=D(i,k) for all primes p, and all 
exponents k needed to reach the preset limit. It then proceeds to factorize consecutive 
values of n. If n is prime then S(n) =n otherwise p,* is replaced by D(i,k) whenever k> 1. The 
largest component in the resulting array is determined and is equal to S(n). Slighly different 
versions of the program has been used depending on the application. Up to n=32000 both 
S(n) and D(i,k) were registed on files. The values of S(n) for n< 4800 were listed with the 
help of a program SN_TAB and the values of S(p;) were listed for p,< 73, k < 75 with the 
help of a program SNP TAB. 


Problem (i) 

(a) No closed expression for S(n) has been found 

(b) No asymptotic expression for S(n) has been found. The behaviour of S(n) for n<32000 
has been graphically displayed using a program SN DISTR. 


Problem (ii) 
S(n)/(n-m), m0 is equivalent to 


m=n+kS(n), k integer .... (1) 
Let us assume that m is a given prime p. From the definition of S(n) it is evident that for 


every n there exists a prime q such that q/n (or n=1q, | integer) and S(n) =jq (i integer). We 
can therefore write (1) in the form 


p=q(1+kj) ... (2) 


To find solutions to (1) when m is a prime pp it is therefore sufficient to chose n as a 
multiple of p which fills the condition 1+kj=1. 


In practice (as far as I have found) this means excluding from n those multiples of m(=p) 
which are divisible by primes larger than p and also cases where n-m has a different parity 
from S(n) as for example (n,m,S(n)) =(5054,19,38) is not a solution while (2527,29,38) is a 


solution. 


When m is not a prime let p be the largest prime such that p/m, i.e. m=rp. Solutions to (1) 
will then be found when n is a multipe of p for which the GCD of n and S(n) is rp. 


The above conditions are sufficient but may not be necessary. Lists of solutions are however 
easily obtained (not included) by looking for solutions to (a-m) mod S(n) = 0. 


Problem (iii) 

A number of solutions to S(x")+S(y*) =S(z") has been obtained and listed for n=3,5,7 and 
11. The program SMAR iii uses only Smarandache function values of the type S(p,*) 
which had first been sorted in ascending order using a program SNP_SORT. 


Problem (iv) 
A program SMAR_iv has been designed to find solutions to the equation S(k*)'=S(n*) but 
no non-trivial solutions were found in the selected search area n< 8000. 


Problem (v) 

In a first attempt values saved on file up to 32000 were in a program SMAR y¥ to calculate 
the sums Z(n)=1+1/S(2)+1/S(3)+ ... +1/S(n) for n= 800, 1600, 2400, ... 32000. These 
sums were used to study the behaviour of Z(n)-T(n) for various functions T(n): 

T(n) =log(S(n)) gave a curve "parallel" to Z(n). 

T(n) =log(largest prime <n) gave a similar result. 

T(n)=14+1/27+1/37+ ... +1/n* gave interesting results. Supplemented with a linear term 
a "nearly straight horizontal” line was obtained. 

To see if this holds for larger values the exercise was repeated for n< 1000000. Computer 
files to store S(n) is now out of question and the generating program SMARAND was 
revised so that the partial sums Z(n) were calculated in the same program. T(n) was 
calculated in a separate program for vaious values of a. For a=0.5 and it was found that 


1+1/S(2)+1/S(3)+ .. +1/S(n) - (1t+1/V2+1/V3+ .. +1/Vn) - (20k-58) 


where k= n/25000 deviates from 0 with at most 10 in the interval 1 to 1000000 (at the 
points of representations in the graph, 1000000 was divided in 40 interval of 25000). 


Problem (vi) 
Not attempted. 


Equipment 

Calculations were done on an dtk-computer with 486/33 Mhz processor. Programs were 
written in Borlands Turbo Basic. Printouts were done on an HP ITP Laser printer. Some 
graphs were done on an HP Paintjet. The run time to calculate S(n) up to 1000000 was 2 
h 50 m. The initial calculation up 32767 took 198 s. 
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*SMARAND, H. Ibstedt, 930320 


‘The Smarandache function S(n) calculated by comparing largest prime and S(P“A). The 
values of S(n) are calculated and registered in a file SN.DAT up to n = 32000. The calculation 
goes further in other applications. 


DEFLNG A-S 
CLS :T=TIMER 

DIM P(168),D(168,75),K(168),L(168) 

OPEN "PA" FOR INPUT AS #1 

FOR |=1 TO 168 :INPUT #1,P(!) NEXT :CLOSE #1 


‘This part of the program calculates S(P(I)~A) and saves the result in the array D(I,A), P(I) is the 
Ith prime number. The routine uses the fact that D(I,A) < =P(I)*A for a downward search for the 
value of D(I,A). This calculation goes beyond what is required to calculate S(n) up to n = 
32000. 


FOR |=1 TO 42 
A=2 :P=P(1) :D(I,1)=P 
WHILE A< 76 
C=0 :N=0 
ke 
C=C+1 
N=N+P 
IF C> =A THEN D(1,A}=N :GOTO LWEND 
PP=P*P 
L1: 
IF N-PP*INT(N/PP) =0 THEN C=C+1 :PP=PP*P :GOTO L1 
IF C> =A THEN D(_A)=N :GOTO LWEND :ELSE L 
LWEND: 
INCR A 
WEND 
NEXT 


‘The array D(I,A) is stored in a file SNP.DAT for future use. 


OPEN “SNP.DAT™ FOR OUTPUT AS #2 
FOR |=1 TO 42 :FOR J=1 TO 75 
PRINT #2,P(1),J,D(1,J) 

NEXT :NEXT :CLOSE #2 


This part of the program calculates S(N). It calls on the subroutine NFACT to express N in 
prime factor form. Factors P(I)*A with A>1 are replaced by D{IA) and placed in array L() 
together with the factors P(I) of multiplicity 1. S(N) is then the largest component of L(). S(N) is 
stored in a file SN.DAT. 


N=1 

OPEN "SN.DAT™ FOR APPEND AS #3 
WRITE #3,1 

WHILE N<32000 

INCR N :print n 

‘Factorize N. 
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GOSUB NFACT 
IF K(0Q)>0 THEN S=P(0) :GOTO LWR 


‘Construct L(). 


FOR |=1 TO 168 :L(1)=0 :NEXT 
C=0 

FOR 1=1TOM 

INCRC 

IF K(I)=1 THEN L(C) =P()) 

IF K()>1 THEN L(C) =D(,K())) 
NEXT 


Find the largest value of L() and hence S(N). 


$=0 
FOR l=1TOC 
IF L()>S THEN S=L() 
NEXT 

LWR: 
WRITE #3,S 
WEND 
CLOSE #3 
T=TIMER-T :PRINT T 
END 


‘Subroutine for factorization of N. 


NFACT: 
FOR !=0 TO 168 :K(!) =0 :NEXT :P(0) =0 
N1=N :M=0 
FOR |=1 TO 168 
LA: 
IF N1-P(I)*INT(N1/P(1)) =0 THEN K(I)=K(I) +1 :M=I:N1=N1/P(l) :GOTO LA 
IF N1=1 THEN |=168 
NEXT 
IF N1>1 THEN P(0)=N1 :K(0) =1 
RETURN 
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"SN_TAB, H. Ibstedt, 930321 


‘This program uses the results stored in the file SN.DAT produced by the program SMARAND 
to tabulate the first 4800 values of the function S(N). 
‘Set [=21 and NB =82 on HPIIP. 


DEFINT I-P,S :DIM S(4800) 
CLS :WIDTH "LPT1:",120 :S(1)=1 :T=TIMER 
OPEN “SN.DAT" FOR INPUT AS #1 

FOR |=1 TO 4800 

INPUT #1,S(1) 

NEXT 


CLOSE #1 
S1S= 1 pelt: SZ SM SBS Ny 
S4$s"] on | Sn) [" :S58=" nm | S(m) [* :86$=" on | S¢n) J" 
S7$=" }—— fit 2 $882 2 $9854 
BUSH" LN BOG EI BBG a a ed 3 
i1=1 :12=75 :P1=1 
LW: 
LPRINT TAB(8) “The Smarandache Function S(n). 
LPRINT TAB(8) S1$; :FOR !=1 TO 6 :LPRINT S2$; :NEXT :LPRINT S3$; 
LPRINT TAB(8) S4$; :FOR |=1 TO 6 :LPRINT S5$; :NEXT :-LPRINT S6$; 
LPRINT TAB(8) S7$; :FOR |=1 TO 6 :LPRINT S8S; :NEXT :LPRINT S9$; 
FOR |=!1 TO 12 
LPRINT TAB(8} “="; 
FOR J=0TO7 
LPRINT USING "#####"}1 + J*75; :LPRINT “™; :LPRINT USING "#####";S(1+J*75); 
LPRINT “oO ; 
NEXT 
NEXT 
LPRINT TAB(8) B1$; :FOR |=1 TO 6 :LPRINT B2$; :NEXT :LPRINT B3$; 
LPRINT TAB(8) “Page"P1 “of 8. 
LPRINT CHR$(12) 
iF P1=1 THEN P1=2 :I1 =601 :12=675 :GOTO LW 
IF Pt =2 THEN P1 =3 :11=1201 :12=1275 :GOTO LW 
IF P1=3 THEN P1 =4 :11=1801 :12=1875 :GOTO LW 
IF P1=4 THEN P1=5 :l1 =2401 :12=2475 :GOTO LW 
IF P1=5 THEN P1=6 :I11 =3001 :12=3075 :GOTO LW 
IF P1=6 THEN P1=7 :11 =3601 :12=3675 :GOTO LW 
IF P1=7 THEN P1=8 :i1=4201 :12=4275 :GOTO LW 
PRINT "END" :END 
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The Smarandache Function S(n). 
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The Smarandache Function S(n). 
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The Smarandache Function S(n). 
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The Smarandache Function $(n). 











The Smarandache Function S(n). 
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The Smarandache Function S(n). 
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The Smarandache Function Sin). 
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The Smarandache Function S(n). 
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"SNP_TAB, H. Ibstedt, 930322 
‘This program tabulates S(n) forn=PJ, 1<P<74, 1<J,76, using data form the file SNP_ASC. 


DEFLNG I-S :DIM KP(1575),KJ(1575),SP(1575) 
CLS: 

WIDTH "LPT1:",120 

OPEN “SNP.DAT™ FOR INPUT AS #1 

FOR l=1 TO 1575 

INPUT #1,KP(I),KJ(I), SP(I) 


NEXT 
CLOSE #1 
s$1$=" end Se es eS a 


S4$="] PCI) | J (DCT, J) J" s85$=" PCT) | os OCI, u)]™ 

a se Oe 

BS ln 2 BOSH et 
a aa caer 

s6$=" PCI) | J [OCI d24" 

2 —— 

fe en FD 

W=1 :12=75 :P1=1 

LW: 

LPRINT TAB(12) "The Smarandache Function S(n) =D(1,J) for powers of primes P()7J.” 
LPRINT TAB(12) S1$; :FOR !=1 TO 3 :LPRINT S2$; :NEXT :LPRINT S38; 
LPRINT TAB(12) S4$; :FOR !=1 TO 3 :-LPRINT S5S; :NEXT :LPRINT S6$; 
LPRINT TAB(12) S7$; :FOR |=1 TO 3 :LPRINT S8$; :NEXT :LPRINT S9$; 
FOR l=11 TO 12 

LPRINT TAB(12) "|"; 

FOR J=0TO4 

LPRINT USING “#####°;KP (1 +J*75); 

LPRINT* |"; :LPRINT USING "#####";KJ(1+J*75); :LPRINT* |"; :LPRINT USING 
“PEALE SP (1+ J5*75); 

LPRINT™ ["; 

NEXT 

NEXT 

LPRINT TAB(12) B1$; :FOR l=1 TO 3 :LPRINT B2$; :NEXT :LPRINT B3$; 
LPRINT TAB(12) "Page"P1 “of 3. 

LPRINT CHR$(12) 

IF P1=1 THEN P1 =2 :11=601 :12=675 :GOTO LW 

IF P1 =2 THEN P1 =3 :11=1201 :12=1275 :GOTO LW 

PRINT "END" :END 
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The Smarandache Function S(n)=0(I,J) for powers of primes PCI) "J. 
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The Smarandache Function S(n)=0(1,J) for powers of primes P(I) “J. 

























































Poy + poop] + powpro] a pow] peop] a pow 
1 3 2 1 rs] 31 1 3 37 1 37 1 41 
2 44 2 58 31 2 $2 37 2 % 2 82 
3 3 3 3 93 37 3 3 
4 4 31 4 124 37 4 4 
5 5 31 5 155 37 5 5 
6 6 31 6 186 37 -) 6 
7 7 31 7 217 37 7 ‘4 
8 8 31 8 2468 37 8 8 
9 9 3 9 279 37 9 9 

10 31 10 310 37 
| 319 31 11 341 37 
12 348 31 12 372 37 
13 377 31 13 403 37 
14 406 31 14 434 37 
15 435 31 15 465 37 
16 bbs 31 16 496 37 
17 493 31 17 527 37 
18 522 31 18 558 37 
19 34 19 589 37 
31 20 620 37 
31 21 651 37 
31 22 582 37 
31 23 713 37 
3 24 Tht 37 
31 3s 775 37 
31 26 806 37 
3 27 837 37 
31 28 868 37 
31 29 899 37 
31 30 930 37 
31 3 961 37 
31 32 961 37 
31 33 992 37 
34 3% | 1023 37 
31 35 | 1054 37 
+3 36 | 1085 37 
31 37 | 1116 37 
38 


















31 1891 37 
31 1922 37 
31 1922 37 
34 1953 37 


1984 37 
2015 37 
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“9 31 2077 37 
70 31 2108 37 
71 31 2139 37 
v2 31 2170 37 
73 31 2201 37 
7% 3 2232 37 
75 31 2263 37 
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The Smarandache Function S(n)=0(1,J) for powers of primes P(I) “J. 
































1 59 61 1 61 67 1 67 71 1 71 1 3B 
2 61 2 122 67 2 134 2 162 2 146 
3 61 3 183 3 201 3 3 219 
4 él 4 244 4 4 4 
5 61 5 305 5 5 5 
& 61 6 366 6 6 6 
7 61 7 427 7 7 7 
8 61 8 488 8 8 8 
9 61 9 549 9 9 9 
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"SN_DISTR, H. Ibstedt, 930322 


‘The values of S(n) for n<32000 are input from the file SN.DAT and the number of values falling 
into each square of a 40 x 40 matrix are counted and displayed in a graph. An interresting 
pattern is formed by large primes while the bottom layer mainly resulting form composite 
numbers requires two lines in the graph. 


‘Set !=1021 and NB=82 on HPIIP. 


DEFLNG A-S 

DIM C(51,51) 

CLS :WIDTH “LPT1:",130 

LPRINT TAB(50) “The Smarandache function S(n)* 

LPRINT TAB(50) een! LPRINT :LPRINT 
LPRINT TAB(35) “Number of values of S(n) in the interval 800y u S(n) < 800(y+ 1). 
LPRINT 

LPRINT TAB(4)” | y~ 

LPRINT TAB(4)" |" 

K=0 

OPEN “SN.DAT* FOR INPUT AS #1 

WHILE K<32000 

INCR K 

INPUT #1,S 

1=K\800+1 :J=S\800+1 

C(i,J) = CU) +1 

WEND 

CLOSE #1 

FOR J=40 TO 2 STEP -1 :LPRINT USING "###";J-1; :LPRINT "|"; 
FOR i=1 TO 40 

IF C(i,J)=0 THEN LPRINT SPC(3); ELSE LPRINT USING "###°;C(1,J); 
NEXT :LPRINT 

NEXT 

LPRINT USING "###";J-1; :LPRINT "|"; 

J=1:FOR !=1 TO 39 STEP 2 

LPRINT USING "###";C(1,J); -LPRINT SPC(3); 

NEXT :LPRINT 

LPRINT USING “###*;J-1; :LPRINT "|"; 

FOR |=2 TO 40 STEP 2 

LPRINT SPC(3); :LPRINT USING "###";C(1,J); 

NEXT :LPRINT 

LPRINT TABC5) "4"; :FOR [=1 TO 120 :LPRINT "-"; :NEXT :LPRINT "> x" 
LPRINT TAB(6); :FOR !=0 TO 39:LPRINT USING "“###°sI; NEXT :LPRINT 
LPRINT :LPRINT TAB(5) “Intervals: 800x< n < 800(x+1). 

LPRINT :LPRINT TAB(5) “SN_DISTR "DATES 

LPRINT CHRS$(12) 

END 
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The Smarandache function S(n) 





Number of values of S(m) in the interval 800y < S(n) < 800(y+1). 


43 
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29 27 32 





45 51 

46 % 

49 47 

47 33 36 29 
Ve] 50 45 29 37 29 21 29 23 4 
106 52 54 33 37 63 2 47 49 38 35 59 3 

112 57 95 67 85 70 87 69 90 69 82 77 79 78 76 
79 §=6©637 «655830 = 555) 5320510 495 BA 
688 610 565 535 521 504 493 478 


35 30 31 
26 25 25 22 21 18 
GO 22 20 31 33 21 27 28 
2 43 49 40 38 53 34 47 
82 73 88 71 79 74 
465 457 
471, 461455 
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Intervals: 800x <s n < 800(x+1). 


SN_DISTR 03-24-1993 
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’SNP_SORT, H. Ibstedt, 930322 


‘This program inputs the Smarandache function S(n) for powers of primes from the file 
SNP.DAT, sorts S(n) in ascending order and writes the result to a file SNP_ASC 


DEFLNG A-S 
CLS 

DIM D(42,75),KP(3150),KJ(3150),SP(3150) 
OPEN "“SNP.DAT™ FOR INPUT AS #1 

FOR |=1 TO 42 :FOR J=1 TO 75 

INCR L 

INPUT #1,K1,K2,K3 

KP(L)=K1 :KJ(L) =K2 :SP(L) =K3 

NEXT :NEXT 

FOR l=1 TO L:FOR J=!+1 TOL-1 

IF SP(I) >SP(J) THEN SWAP SP(1),SP(J) :SWAP KJ(l),KJ(J) :SWAP KP(I),KP(J) 
NEXT :NEXT 

OPEN “SNP_ASC* FOR OUTPUT AS #2 
FOR l=1 TOL 

PRINT #2,KP(I),KJ(1),SP(1) 

NEXT 

CLOSE #2 

PRINT "END" :END 
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"SMAR _ii, H. Ibstedt, 930322 


This program searches for solutions to S(x*n) +s(y~n) = S(z*n). Two parameters are set in the 
program: n (=NM) and the largest value of S(z*n) (=NS). x, y and z are restricted to powers of 
prime numbers. The input to the program is provided by the file SNP_ASC, which contains 
S(p*n) sorted in ascending order. 


DEFLNG A-P :DEFDBL X,Y,Z 
CLS :WIDTH "LPT1:",120 

DIM D(42,75),K1(3150),K2(3150),K3(3150) 

NM=3 :NS=120 :L=0 

OPEN "SNP_ASC* FOR INPUT AS #1 

WHILE NOT EOF(1) 

INCR L :INPUT #1,K1(L),K2(L),K3(L) 

WEND :CLOSE #1 :COUNT=0 

LPRINT TAB(16) "Solutions to $(x*n) + S(y“n) = S(z*n) forn ="NM 
LPRINT TABC16) :LPRINT 

fn ee ERE meen REE eee 
—_—" 

LPRINT TABC16) :LPRINT “f PIT | P2 | P3 | x | y | z 
J Sex'n) | Sty°n) | Séz?n) J* 

LPRINT TAB(16) :LPRINT 
‘HH HH Ht 
FOR |=NS TO 4 STEP -1 :PRINT | 

IF K2(7) MOD NM <> 0 THEN il 

FOR J=I-1 TO 3 STEP -1 

IF K2(J) MOD NM <> 0 THEN JJJ 

FOR K=J-1 TO 2 STEP -1 

IF COUNT =75 THEN GOTO KKK 

iF K2(K) MOD NM <> 0 THEN KKK 

IF K1() =K1(J) OR K1(1)=K1(K) OR K1(J) =K1(K) THEN KKK 

IF K3(1) =K3(J) + K3(K) THEN 

INCR COUNT 

X=K1(K)*(K2(K)/NM) :¥ =K1(J)*(K2(J)/NM) :Z=K1 (1)*(K2(1) /NM) 
LPRINT TAB(16) ]°; :LPRINT USING "######";K1(K); 

LPRINT "|"; :LPRINT USING "######°;K1 (J); 

LPRINT "|°; :LPRINT USING "######°K1 (1); 

LPRINT "}°; :LPRINT USING "#######4#4#¢###4#":X: 

LPRINT "|"; :LPRINT USING "############EHCY: 

LPRINT "|"; :LPRINT USING "############4#7:Z. 

LPRINT "|"; :LPRINT USING "#######°;K3(K); 

LPRINT °|"; :LPRINT USING "####4##€#";K3(J); 

LPRINT "|"; :LPRINT USING "#######°;K3(1); 

LPRINT "]*° 

END IF 

KKK: 

NEXT 

JJ: 

NEXT 

UH: 

NEXT 

LPRINT TABC16) :LPRINT 
ee Ce eS ee ee Se ee 
—!" 

LPRINT CHR$(12) 

PRINT “COUNT ="COUNT 

PRINT “END* :END 


Solutions to S{x7n) + Syn) = S(2z°n) for n= 
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2 
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7 
2 
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5 
5 
7 
2 
5 
5 
3 
2 
3 
5 
3 
5 
7 
3 
2 
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Sotutions to S(x°n) + Syn) * S(z°n) for n = 5 
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7 
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5 
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3 
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5 
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1594323 
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64 
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25 
55 
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Solutions to S(x°n) * SCy"n) = S(z°n) for n= 7 
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Solutions to S(x°n) + SCyn) = S(z°A) for n = 


CC 


pellet 
68921 
68921 
68921 
3721 
3721 
3721 
3721 
923521 
3481 
3481 
707281 
707281 
50653 
50653 
50653 
50653 
6436343 
6436343 
2809 
2809 
26137569 
26137569 
24137569 


161051 
1369 
1771561 
61 
53 
161051 
4913 
43 
43 
12167 
1771561 
1771561 
61 
61 
37 
37 
961 
28561 
841 
529 
371293 
28561 
43 
41 
59 
6859 
43 
16641 
41 
289 
529 
47 





"SMAR_iv, H. ibstedt, 930323 


‘This program searches for solutions to the equation S(k*n)*i = S(n“k). The search is limited 
to the first 8000 values of S(n) loaded from the file SN.DAT. No non-trivial solutions were 
found. 


DEFLNG I-S :DEFDBL X,Y,Z 
CLS :WIDTH "LPT1:",120 
DIM S(8000) 
L=0 :NS=90 :KS=90 
OPEN "SN.DAT" FOR INPUT AS #1 
FOR L=1 TO 8000 
INPUT #1,S(L) 
NEXT 
CLOSE #1 
PRINT S(8000) 
LPRINT TAB(12) “Search for solutions to S(k*n)“i = S(n*k). * 
LPRINT TABC12) SLPRINT pene repent nnn 
LPRINT TAB(12) :LPRINT "f k [| mn | Stk*mfo nm [| k | Stk] i |" 
LPRINT TABC12Z) cLPRINT ' fpsnnmenniseennnn een tne heed 1 
FOR K=2 TO KS :PRINT K 
FOR N=2 TONS 
iF K°N>8000 OR N*K>8000 THEN N=NS :GOTO L2 
C=S(K°N) :D=S(N*K) 
IF C>D THEN SWAP C,D 
E=1 :!=0 
L1: 
E=E*C :i=i+1 
iF E>D THEN L2 
IF E=D AND K<>N THEN GOSUB LW 
IF E<D THEN L1 
L2: 
NEXT 
NEXT 
LPRINT TABC 12) SLPRINT "lemme haem 
LPRINT CHR$(12) 
PRINT “END” :END 
LW: 
_ LPRINT TAB(12) "["; :LPRINT USING "####4##4°K; 
LPRINT "|"; :LPRINT USING “###4#4##°;N; 
LPRINT "|"; :LPRINT USING "######*;S(KAN); 
LPRINT "|"; :LPRINT USING “####4##°,N; 
LPRINT "|"; :LPRINT USING “#####4#°:K; 
LPRINT "|"; :LPRINT USING "######°;/S(N“K); 
LPRINT "|"; :LPRINT USING "#######4')1; 
LPRINT "|" 
RETURN 


Search for solutions to S¢(k°n)"i = S(n°k). 
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"SMAR_y, H. Ibstedt, 930323 


‘This program uses the first 32000 values of S(n) input from the file SN.DAT to calculate 
1+ 1/S(2)+ 1/S(3) + ... + 1/S(n) for n=800, 1600, 2400, ... 32000. The sums are used to study 
the behaviour of 1+ 1/S(2) + 1/S(3) + ... 1/S(n)-T(n). The following cases are examined T(n) =0, 
T(n) =Log(S(n)), T(n) equal to the logarithm of the largest prime less than n and finally 
T(n) =(rK-V) + 1/2°at+1/Sa+ ... +1/n°a, where a=.5164, r=0.96 and V=27. K=n/800. 


‘Set [=22 and NB=62 on HPIIP 
‘Number of intervals = NN (40), Size of interval NI (800). 


DEFDBL A-Z :DIM LAM(40),SAM(40),ZAM(40), HAM(40) 

NN =40 :NI=800 

CLS :WIDTH “LPT1:",120 

Z=1 :K=1:21=0 

LPRINT :LPRINT 

LPRINT TAB(12); “Behaviour of 1+1/S(2)+1/S(3)+ ... 1/S{n)-T(n). K=n\800." 
LPRINT TAB(12); “In the column LAM(K) T(n)=0, in SAM(K) T(n) =Log(S(n)),in* 
LPRINT TAB(12); “ZAM(K) T(n) is the logarithm of the largest prime less" 
LPRINT TAB(12); “than n and in HAM(K) T(n) =1+1/2%.52+1/37.52+ ... 1/n%.52." -LPRINT 
LPRINT TABC 12) 7 1 ny yt 
LPRINT TAB(12); "| K | LAMCK) |  SAMCK) |  ZAMCK) |  HAMCK)  [" 
PRINT TAB C12) 5 1 pen en es 
OPEN “SN.DAT™ FOR INPUT AS #1 

WHILE K<41 

G=1 

WHILE | <K*Ni 

INCR | 

INPUT #1,S 

IF S>G THEN G=S 

Z=Z+1/S 

Z1=Z1+1/17.5164 

WEND 

LAM(K) =Z 

HAM(K) =Z-Z1-.96*K+V 

SAM(K) =Z-LOG(S) 

ZAM (K) = Z-LOG(G) 

LPRINT TAB(12); :LPRINT ‘}"; :LPRINT USING "#########4#4':K; 

LPRINT "|"; :LPRINT USING "#####4###4##";LAM(K); 

LPRINT "|"; :LPRINT USING "##########";SAM(K); 

LPRINT "|"; :LPRINT USING "#######4###";ZAM(K); 

LPRINT "|"; :LPRINT USING "##########°7HAM(K); 

LPRINT "4" 

INCR K 

WEND 

CLOSE #1 

LPRINT TAB(12); a 
LPRINT :LPRINT TAB(12) “SMAR_v" 

LPRINT CHR$(12) 

PRINT “END" :END 
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PCW, Numbers Count: Problem {v), — ae 
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1+1/S(2)+ ... +1/S(n)-log(S(n)) 
n = 800, 1600, 2400, ... 32000 


Samma (S) 
g 


8 


S 


1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31.33 35 37 39 
2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 


= n/800 


Behaviour of 1+1/8(2)+1/S(3)+ ... 1/S(n)-T(n). K=n\800. 

In the column LAM(K) T(n)=0, in SAM(K) T(n)=Log(S(n)),in 
ZAM(K) T(n) is the logarithm of the largest prime less 

than n and in HAM(K) T(n)=(r*K-V)+1/2°at+1/3°a+ ... Ifn°a, 
where r=0.96, V=27 and a=0.5164 have been chosen to fit HAM(K) 
as closely as possible to 0. 


























K LAM(K) SAM(K) ZAM (K) HAM (K) 
1 37 35 31 13 
2 55 53 48 9 
3 70 67 62 6 
4 82 80 | 74 5 
5 94 91 85 3 
6 104 102 96 2 
7 114 112 105 1 
8 123 121 115 i 
9 132 130 123 oO | 

10 141 138 132 -0 
il 149 146 140 -0 
12 157 154 148 -1 
13 164 162 155 -1 
14 172 169 162 -1 
15 179 176 170 -1 
16 186 183 176 -1 
ig 193 190 183 -1 
18 200 197 190 -1 
19 206 203 196 -1 | 
20 212 210 203 -1 
21 219 216 209 -1 
22 225 223 215 -1 
23 231 228 221 -1 
24 237 235 227 -1 
25 243 240 233 -1 
26 249 246 239 -1 
27 254 252 244 -1 
28 260 258 250 -1 
29 265 262 255 -1 
30 271 268 261 -1 
31 276 273 266 -1 
32 282 279 272 -1 
33 287 285 277 ~0 
34 292 289 282 -0 
35 297 295 287 -0 
36 303 300 292 -0 
37 307 304 297 -0 
38 313 310 302 fe) 
39 317 315 307 @) 
40 322 320 312 fe) 


SMAR_v 


The Smarandache Function S(n) 


Graphical representation of 1+ 1/S(2)+1/S(3) + ... + 1/S(n) - T(n) for n< 32000. 


Graph |: T(n) =0. 

Graph Il: T(n) =log(S(n)). 

Graph III: T(n) =log(largest prime<n). 

Graph IV: T(n) = (rK-V) + 1/2%a+1/3%a+ ... +1/n“a, where the parameters r=0.96, V=27 and 


a=0.5164 were chosen to fit the graph as closely as possible to a horizontal line. 
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*SMARAND1, H. Ibstedt, 930327 


The Smarandache function S(n) calculated by comparing largest prime and S(P*A). The 
upper limit for the calculation is n = 1000000. The results are used to calculate 
1+ 1/S(2)+ 1/S(3) + ... +1/S(n) registering partial sums for n = 25000, 50000, 75000, ... 
1000000. 


DEFLNG A-S 

CLS :T=TIMER 

DIM P(168),D(168,20),K(168),L(168) 

OPEN “PA" FOR INPUT AS #1 

FOR |=1 TO 168 :INPUT #1,P(I) NEXT :CLOSE #1 


‘This part of the program calculates S(P(I)“A) and saves the result in the array D(A), P(I) is the 
Ith prime number. The routine uses the fact that D(A) < =P(I)*A in the search for the value of 
D{LA). This calculation goes as far as is required to calculate S(n) up to n = 1000000. 


FOR !=1 TO 168 
A=2 :P=P(l) :D(I,1)=P 
WHILE A<21 
C=0 :N=0 
L: 
INCR C 
INCR N,P 
IF C> =A THEN D(1,A)=N :GOTO LWEND 
PP =P*P 
Lt: 
IF N-PP*INT(N/PP) =0 THEN INCR C :PP =PP*P :GOTO L1 
IF C> =A THEN D(I,A) =N :GOTO LWEND :ELSE L 
LWEND.: 
INCRA 
WEND 
NEXT 


‘This part of the program calculates S(N) and the sum of reciprocals. !t calls on the subroutine 
NFACT to express N in prime factor form. Factors P(I)“A with A> 1 are replaced by D{I,A) and 
placed in array L() together with the factors P(I) of multiplicity 1. S(N) is then the largest 
component of LQ). S(N) is stored in a file SN.DAT. 


N=1:Z=1 :D=0:ZC=0 

OPEN “SAM.DAT™ FOR APPEND AS #3 
WHILE N< 1000001 

if inkey$< >"" then print “end" :end 
INCR N :INCR D :PRINT N 

‘Factorize N. 

GOSUB NFACT 

IF K(0)>0 THEN S=P(0) :GOTO LWR 
‘Construct LQ). 

FOR |=1 TO 168 :L(l) =0 :NEXT 

C=0 

FOR |=1 TOM 

INCR C 


IF K(l)=1 THEN L(C) =P() 
IF K()>1 THEN L(C)=D(1,K(I)) 
NEXT 


Find the largest value of LQ and hence S(N). 


$=0 
FOR !=1TOC 
IF L(l)>S THEN S=L(I) 
NEXT 

LWR: 
Z=Z+1/S 
IF D=25000 THEN ZC =ZC+Z :WRITE #3,Z,ZC :Z=0 :D=0 
WEND 
CLOSE #3 
T=TIMER-T :PRINT T 
END 


‘Subroutine for factorization of N. (Improved from SMARAND to avoid large primes) 


NFACT: 
FOR I=0 TO 168 :K(1)=0 :NEXT :P(0)=0 
N1=N :1=0 :M=0 
FOR |=1 TO 168 
LA: 
IF N1-P(1)*INT(N1/P(1)) =0 THEN K(1) =K(I) +1 :M=I :N1=N1/P(l) :GOTO LA 
iF N1=1 THEN 1=168 
NEXT 
IF N1>1 THEN P({0)=N1 :K(O}=1 
RETURN 
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Problem v: Summary of obtained data. 


Graph | Graph Il 
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The Smarandache Function S(n) 
Comparison with the sum of 1/n* 0.5 


for n up to 1000000 





13 5 7 9 1113151719 21 23 25 27 29 31 33 35 37 39 K 
2 4 6 8 1012 1416 18 20 22 24 26 28 30 32 34 36 38 40 
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